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ABSTRACT 


Electron velocity distributions are deduced along axes parallel and 

perpendicular to the magnetic field in a pulsed, linear Penning discharge 

in hydrogen by means of a laser Thomson scattering experiment . Results 

obtained are numerical averages of many individual measurements made at 

specific spacetime points in the plasma evolution, and take advantage of 

the high level of reproducibility of density and velocity parameters 

possible in this type of discharge. Observation solid angles are made 

quite small to maximize resolution between the two distributions. Because 

of the high resolution in k-space and the relatively low maximum electron 
1 ^ -3 

density 2 X 10 cm , special techniques were required to obtain measurable 
scattering signals. These techniques are discussed and experimental A 
results are presented. 

Our Thomson scattering measurements show identical behavior of 
the longitudinal and azimuthal distributions, within experimental error. 

For magnetic fields below the critical value for onset of strong insta- 
bility in the Penning discharge, bur results show impressively good fit 
to a calculated gaussian distribution, for electron energies up to 6 
times the thermal energy. However, measurements above this critical 
field show a distinct excess of high energy electrons which cannot be 
accounted for by a single gaussian, within experimental error. Equipment 
scattered light levels are the lowest ever attained in an experiment of 
this kind. 
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1 . INTRODUCTION 


1 .1 General 

Since the availability of ruby lasers made Thomson scattering 

1 2 

practicable in laboratory plasmas ’ this technique has attracted 
intense interest as a plasma diagnostic tool. The prime reason for 
this interest is that Thomson scattering makes directly available the 
ion and electron velocity distributions and thus in principle does not 
require thermal equilibrium for either species. Furthermore it can 
be used anywhere a beam of light^ can be sent, provided the electron 
density is sufficient to produce a detectable scattered signal. In 
particular, this method is well suited for regions of density and 
temperature inaccessible to standard techniques, such as occur in thermo 
nuclear plasmas. 

In this work, we are reporting results from an experiment in 
which the electron velocity distributions parallel and perpendicular 
to the axis of a linear Penning discharge plasma are measured. Our 
purpose is to determine whether the overall distribution is anisotropic. 
For the first time, these distributions are measured separately and un- 
equivocally at a particular spacetime coordinate, with convincing in- 
dependence between the two measurements. The Penning discharge is 
used because it can be conveniently driven from quiescent operation 
to strong drift-type instabilities at sufficiently high density to 
make scattering measurements feasible. We operate the machine in a 
pulsed mode so that the dynamic history of velocity and density para- 
meters may be constructed, as a function of the applied magnetic field. 
The discharge is operated in hydrogen. 

This experiment is a member of the emerging second generation of 
Thomson scattering experiments, in which this technique is used as a 
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well understood tool and where the surprises are expected from the plasma 


rather than from the physics governing scattering. 

For the benefit of the reader who has not used this technique , it 
should be mentioned that such experiments invariably involve a large 
initial investment of time (in our case approximately * 10 man hours) , 
considerable expense and patience. Nevertheless, it is a powerful, pre- 
cise and widely applicable method which can yield measurements of plasma 
parameters that can be obtained in no other way. 

1 .2 Brief Introduction to Thomson Scattering and Its Applications 

Electromagnetic wave scattering by density fluctuations in a plasma 
has been properly understood for just over a decade. Prior to the 

*3 

pioneering experiment of Bowles^ in 1958; the richness of information 
contained in the scattered spectrum was not realized. This experiment 
was the stimulus for intensive theoretical work in the early 1960 ' s^ ^ 
which resulted in a complete scattering theory based on a modern under- 
standing of plasma fluctuations. 

4 

Gordon suggested the Bowles experiment in a paper which showed 
the feasibility of using radar backscatter measurements for the remote 
determination of electron density and temperature at any level in the 
earth’s ionosphere. It was originally thought that the backscattered 
microwave spectrum would have a thermal doppler width characteristic of 

V 

the electron temperature at the chosen level . The expected spread 
of frequencies in the scattered signal was of the order of 20 Khz. at 
a radar frequency of 40 Mhz. However, Bowles' results showed a spread 
many times smaller than could be explained by the lowest available 
estimates of electron temperatures at the observation height. 
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Bowles himself realized that his results were explained if the 
observed scattering was produced by just the electron density fluc- 
tuations required to neutralize the slow density fluctuations of the 
more massive ions in the plasma. Of course, the electron spectrum is 
not always absent, and with the benefit of modem scattering theory we 
would say that a particular combination of frequency, scattering angle 
and electron density and temperature permitted Bowles to observe only 
the "ion feature” of his plasma (see Chapter 2), In contrast, our 
experiment is concerned only with the "electron feature", where the 
scattered signal is doppler shifted according to the electron velocity 
distribution. 

To investigate a plasma in the laboratory, it is necessary to scale 
down the illumination wavelength to the optical region, for two reasons. 
First, laboratory plasmas tend to have Debye lengths of the order of 
several microns at most and, as will be seen in Chapter 2, the illu- 
mination wavelength must be much less than times the Debye length 
in order to observe electron features. Second, any experiment designed 
to measure scattering off thermal fluctuations in a plasma requires 
quite large photon densities in the illumination beam just to get one 
or two photons into the detector. It happens that lasers are the only 
adequately bright primary beam sources of appropriate wavelength 
available today. Visible or near infrared lasers are indicated because 
detectors are more sensitive in this region. In plasmas with very 
large Debye lengths where far infrared or microwave sources would be 
permitted in principle, the diffraction limited spacial resolution is 
insufficient for the small plasmas found in laboratories. 

It is true that powerful x-ray sources are also available. Un- 
fortunately, these are not sufficiently bright or monochromatic, and 
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individual x-ray photons are so energetic that, for fixed illumination 
energy, fewer photons enter the experiment. At the detector, this 
is a statistical disadvantage. 

In the next chapter, we discuss the theory of Thomson scattering, 
before continuing the history of scattering experiments in Chapter 3- 



2. THOMSON SCATTERING - THEORY 


2 .1 Definition of Scattering Types 

Photon " scattering' 1 occurs when a beam of electromagnetic radiation 
penetrates a material medium. The electric field of the illuminating 
wave produces microscopic electric currents which are hen the source of 
secondary or ” scattered” photon radiation. The latter will generally 
travel in different directions and have different frequencies than the 
illuminating wave. 

’’Thomson" and !, Compton” scattering are subclasses in which the motion 

17 

of free charged particles causes the scattering. Thomson scattering 
is the long-wavelength, classical limit of electric monopole scattering, 
and applies when the photon energy is much less than or approximately 
equal to the particle rest energy. 

In contrast, when the scattering is due to the motion of bound elec- 
1 8 

trons, Mie-Debye scattering is obtained. Rayleigh scattering is the 
long-wavelength limit of this process, for frequencies well below resonance 
for the electrons in a molecule or ion system. 


2 .2 Scattering by Single Nonrelativist ic Electrons 

According to the Larmor formula, a free electron being accelerated 
by the electric field of a plane, linearly polarized, monochromatic 
electromagnetic wave E - E^e^^i X gives rise to a reradiated field 

of instantaneous magnitude 


e 2 E 


r ! 

me R 


E r 

° sinY = ~ ■ sinY 


( 2 - 1 ) 


** 


P 

For scattering from electrons, with m e c^ = 511 Kev, the Thomson approxi- 
mation is valid up through the x-ray region. Thomson* s work was in fact 
done to describe the scattering of x-rays by electrons. 

Gaussian units are used throughout this work, except as noted. 
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in direction Y relative to and distance R, so that average 

power into unit solid angle around this direction is 




r 2 sin 2 Y . 
o 


( 2 - 2 ) 


2 

If this power is normalized to unit incident intensity cE q /8?t, the 
’’total Thomson differential scattering cross-section” is obtained, 



2 • 2 ™ 
r sm Y 
o 


(2-3) 


where r^ is the square of the classical electron radius and equal 

«2 6 2 

to 7.94 X 10 cm . We see that Thomson scattered radiation, like any 
dipole radiation, is isotropic in the plane normal to and linearly 

polarized in the plane containing E^ and R. The power of Thomson 
scattering as a diagnostic method lies in the extreme simplicity of 
this relationship between the incident and scattered fields. It is vio- 
lated only when the illuminating field is intense enough to cause rela- 
tivistic motions of the electron, but to avoid this it is sufficient to 
have 

EX « jt(m c 2 /e) = ttW volts , (2-4) 

o o v e ' o 

where W q is the electron rest energy. This is an easy limit, not ex- 
ceeded e.g. for 0.7u radiation until E = 2 X 10^ v/cm or I = 5 X 10 ^ 

1 ^ o o 

2 * 

watts/cm . 


This criterion is mathematically the same as Bernstein’s requirement 
that the electron’s excursion during illumination be less than a wave- 
length of the incident wave. 
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2 . 3 Scattering by Cloud of Independent Electrons 


For a tenuous cloud of randomly spaced electrons (but no ions as 
-1/3 

yet) such that n^ J there will hardly ever be more than one 

electron in a cubic wavelength cell, and their vibrations in the illuminating 
field will be uncorrelated . Then the scattered photon number from a volume 
V of such electrons into solid angle dfi will be 


N /d<j \ 

= n vf-^J dn 

N_ e \dQ / 


( 2 - 5 ) 


with N the number of incident photons per unit area. This is just the 
o 

sum of the individual radiated intensities. 

If the cloud is still composed of independently moving randomly 

-1/3 

spaced electrons, but is more dense such that n^ « X q , we may con- 
sider the total scattering for the sake of discussion to be the sum of 
contributions from adjacent cells in the volume V, each cell having 
edge length X q / 2. If the density were absolutely uniform, at any given 
instant there would be no net scattered electric field because of exact 
phase cancellation between adjacent cell scattered fields. If we allow 

random (thermal) density fluctuations, however, say a number (N. + 6N.) 

0 J 

of electrons in the jth cell, this cell’s contribution to the scattered 
field will be proportional to 6Nje i( ^j and the ensemble averaged total 
scattered intensity to 


("s.) * ^“jl 2 ) 

= <5 (8 V V <£ ji* el( * J '^S 6, 0 

= <5 (6N J )2 > = ? N j = n e y . 


( 2 - 6 ) 


This is because the assumption of independence requires the density 
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fluctuations in adjacent cells to be independent and thus <^ 6 N 6 N^ = 

^ 6 N.^ 6 N k ^ = 0. Therefore (2-5) describes scattering from a cloud of 

3 * r 

independent electrons whether or not n \ > 1 . 

This discussion shows the importance of density fluctuations to 
scattering and predicts that the calculation of scattered power will 
be complicated by the absence of independence, as for example when coherent 
density waves are present. In fact, in the limit of perfect dependence, 
the term we ignored in ( 2 - 6 ) becomes n times larger than the one we 
kept, where n is the number of cells in V, according to the present 
simple analysis. In actuality, the coherent cross-section is arbitrarily 
large . 

2 . 4 Scattering by Plasma: Heuristic Results 

The introduction of an equal number of positive ions into the elec- 
tron cloud produces a much more complicated situation in which the elec- 
trons are less independent and scatter less efficiently ’ The 

primary effects are on the frequency distribution of the scattered light 
which we will consider in the next section. It should be mentioned that 
equation ( 2 - 1 ) shows that scattered power varies inversely as the mass 

squared of the scattering particle, so that direct scattering by the ions 

-7 

of a hydrogen plasma should have an intensity only 3 X 10 of that due 
to the electrons, and less for heavier ions. But the ions do have a 
strong effect on the scattered spectrum through the electron velocity 
distribution . 

We define the important parameter OC as the ratio of the "momentum 

2 1/2 

exchange length" 1/k to the electron Debye length = (KT /4rtn e ) 


If, however, the electrons are neutralized by even a smeared-out 
non-discrete positive charge background, more exact calculation are 


necessary 


12 


8 



* 

where k = Ikl = |k -k. I = 2k. sin and “hk is the amount of momentum 
11 si' 1 2 

exchanged when an incident photon with momentum tik^ is deflected through 
a scattering angle 0 to the direction k g of the scattering observer. 
Then 


a = 


IkDj 


(2-7) 


15 


As shown by Buneman , the overall effect of the ions on the total 
scattering cross-section can be expressed by a modulation function S(a) > 
so that equation (2-3) for electron cloud scattering is replaced by: 



S(a)r o 2 sin 2 Y 


(2-8) 


and 



(2-9) 


where Z is the ionic charge number. This result is more or less inde- 
pendent of the exact nature of the electron and ion velocity distribution 

functions so long as an approximate temperature exists. The same result 

12 

was obtained also by Salpeter . The main features of this modulation 
are that it tends to ( , , ^ 1 in the limit a » 1 and 1.0 in the 

yl+ZTg/Ti J 

limit ot « 1 (see Figure l) . The case a « 1 will be seen to be equi- 
valent to the case of free independent electron scatterers in what is to 
follow. 


For nonrelativistic plasmas, the doppler shift of scattered light, 
being proportional to v/c, will be only a few percent. Thus the 
change in photon momentum magnitude |k|/k i is ignorable in most 
cases . 
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Figure 1 

Total Thomson Cross-Section Modulation Function 


Unfortunately a more accurate calculation of S(a) is necessary 

for large OL and large temperature ratios T^/T^ and especially when 

large electron-ion drift velocities are present as well, due to the ion- 

a> fZKT~ 


acoustic resonance at frequencies 


corresponding to the ion 


sound speed. This resonance can contribute arbitrarily large scattering 

lg , x 

analogous to critical opalescence in solids , whereas (2-9) predicts 

zero scattering under these conditions. Also, we need a calculation that 

shows the frequency dependence of Thomson scattering as well as the k- 

dependence we have been considering up to now. 

* 

At the beginning of this section, we provided an initial motivation 
for the consistent appearance of the momentum difference vector hk, 


rather than the incident or scattered beam vectors hk. or hk in the 

1 s 
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scattering formalism, by mentioning momentum conservation. In fact, a 
moment’s heuristic thought will make the qualitative predictions of the 
scattering theory quite plausible. If we observe (in the scattering 
plane) a photon scattered into a small solid angle around direction 0 
to the illumination direction, we are seeing the result of momentum 
transfer from a plasma electron to the incident photon in the direction 
k ( see Figure 2) . 



Ink = ft (k s - kj) 


Figure 2 

Illustrating Momentum Transfer 

The amount of the electron’s total momentum m v which is transferred 

e 

in a particular event depends only on its momentum vector component along 
k just prior to illumination. If the detector has frequency discrimi- 
nation ability, the amount of momentum transfer can be measured quite 
easily from the scattered photon’s doppler shift 0) = a) g -a^ = k • v^ 
when the electrons are nonrelativistic . A calculated modulation func- 
tion S(k,o)) (also called the dynamic form factor) determines the extent 
to which density fluctuations attenuate or enhance the scattered photon 
flux in a particular plasma. However, it will be shown that this func- 
tion is simple in the absence of large scale coherent density waves. 
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If we momentarily ignore the effect of this modulation (independent) 
electron scatterers, section 2.3), the measured scattered photon inten- 
sity into solid angle dQ around a fixed direction 0 will be seen 
to have a frequency distribution 

dcr 


N dQ ^ek 


(.) 


dQ 


( 2 - 10 ) 


which is proportional to the one-dimensional electron velocity distribu 
tion function f 


ekl 


along the direction k. In (2-10), this function 
is normalized so that its integral is unity. The primary effect of 
varying 0 will in this case be to change k and thus the velocity 
scale factor of the observed distribution. The only additional limita- 
tion that arises in this discussion is the implicit assumption from which 
the power of Thomson scattering techniques derives, that the amount of 

an electron’s total momentum m v transferred to a scattered photon is 

e 

very small. This will be true for any 0 if h/\ « m v or, which is 

o e 


the same thing if 


hv 

— » - 

c ^ 2 

m c 
e 


( 2 - 11 ) 


This is an easy limit also, since for 0.7p radiation it amounts to 

avoiding measurements on electrons with temperatures lower than 0.05°K. 

It is extremely unlikely that an electron will emit more than one scattered 

photon in a particular scattering event, since the probability of a 

6 7 x io~ 2 5 

particular electron scattering one photon is only — — — , where 


A is the cross-section area of the illuminating beam in cm . 


2 .5 Scattering by Plasma: Direct Calculation 

20 

Since Akhiezer et al . first indicated the importance of thermal 
density fluctuations in the calculation of plasma Thomson scattering, 
several authors have contributed to an exact calculation and useful 
approximations which we outline here^ )< ^ >y 21 . 
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We consider a volume V of plasma containing N = n^V electrons 

with V » D . From an origin within V, we will note the observer 
e 

position by R and the position of the jth electron by x.. Under 

__ j 

the influence of an illuminating wave E = E e^ k i X the 

o 

scattered electric field due to acceleration of the jth electron is: 


E (R,t) [R x (R X i)] . (2-12) 


sj 


2 

me R 


R-x.' 1 


, 1 j 1 — / — / 

In the above, the retarded time tj = t — , xj = xj(tj ), 

1 


denotes a unit vector and we have made the approximation 


R-xh 


1 

R 


in the field amplitude bnly for an observer located well away from the 
scattering volume. Now we approximate the retarded time by 


R • x (t-R/c) 

t. 7 « t - S + J 

J c c 


(2-13) 


so that the phase factor becomes -approximately 


» exp i{k 


. ■ *.(• - ?) 



r / \ 

\ 


t . * + iiiliil) 


/ ° 

C c 

( 


- 

) 


exp i{k • x 


L R \ l R 

j \ C / o c_ 


(2-U) 


where we have introduced the difference vector k = k -k 

i s 

Therefore, 


CD R 

k. - — 
1 c 


E 


sj 


( B ’ t + i)“ 


it o 


(2-15) 


Now, the Thomson cross-section, which is just the total scattered 

A 

intensity per unit solid angle in direction R normalized to unit 
incident intensity, can be written (since |r X R X E^j = sinY) 
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(2-16) 


d °T I .( E ' t + f) _ " 2 <|E 3 ,(H,t)| 2 > 



o 



e i(k-Xj(t)) 



In the above, y indicates an ensemble average over all states of the 

plasma (assumed stationary), x^ = x^(t) and the first Born approximation 

(|e ] « |e |) is assumed to be valid. We note that if we define the 
s o 

density of point electrons in the plasma by 


_ N _ _ 

n e ( x ) = 5( x “ x j) (2-17) 

W 

at time t, the space Fourier transform of this quantity is 


So, 


n e (k) . Y) *' lk ' X J (t) . (2-18) 

j^l 

a£- r o 2 ,ln 2 , <l n .< E) l 2 > • (2 " 19> 


If we now also require the frequency dependence of scattered radiation, 
the same procedure yields 


da 

T 

dOdoo 


= r ^sin*TT S(k,ao) 


( 2 - 20 ) 


S(k,co) = |n e (k,a ))| 2 


( 2 - 21 ) 


lb 



%(“<“>> - ““hi 1 


dt e 


icut 


-ik*x (t) 
e j v ' 


( 2 - 22 ) 


From this derivation and the Parseval theorem, we also have 


S(k,u>) = J* dx J dt e i ^ k X ,a)t ^C(x,t.) 


(2-23) 


V t 


where C(r,t) = (n(x,t , )n(x+r,t , +t))- / is the "time dependent pair 

correlation function” of density fluctuations in the plasma. 

Equation (2-20) forms the basis for all scattering calculations. 

The dynamic form factor S(k,co) contains all information about the plasma 
density fluctuation spectrum, and its calculation is only symbolically » 
possible in the general case. Fortunately, useful approximations have 
been developed . 

19 

In a typical approach , the linearized Vlasov-Poisson equations are 
used within the random phase approximation to compute the response of a 
plasma to a test charge. This response is then averaged over all particles 
in the plasma, considering each to be a test charge to obtain 


S(k,o>) = 


1 + 


J *itX 1 (kj<»>) 


«(k,aO 


.(s> 


f e m + z 


^^(kjoo) 

«(k,u) 


'i(f) 


(2-24) 


where 


J d ternr)j 


(2-25) 


6 | 

(Z=l for electrons) are the free electron and ion polarizabilities, and 


m k 
e 
i 


c(k,o>) = 1 + 4jr[X e (k,co) + ^ (*><■>)] j 


(2-26) 


is the plasma dielectric function. 
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In the Salpeter approximation, it is assumed that the velocity dis- 


tributions are Maxwellian with electrons and ions allowed to have different 
temperatures. In this case, Salpeter finds essentially 

v2 


with 


S(k,cn) 


1 + 




2 -(a/u> e )‘ 

e 


'/it 


OD 


+z (5) 


ZT / 2 


1 + 


'(*)•© 


2 -(ay'oo i )‘ 

e 


\/it cu. 


(2-27) 


and 


W(?) = 


g2 

i 

OJ 

1 

1 - 2%e~^ 

1 e P dp - i\£x 5 e i 

< 

o ) 


(2-28) 


(D 


e >i _ 

k "\l 


2KT 


e, l 


m 


e, i 


(2-29) 


is the electron or ion thermal velocity. For OC « 1, the Maxwellian 
dynamic form factor reduces to 


lim 

a-o 


S(k,cjo) 


-(cr/a ) e ) 2 

e 


(2-30) 


the so-called ’’electron feature”. In the opposite limit OC » 1, the 
electron contribution condenses into 2 narrow Lorentzian spikes symmetri- 
cally displaced from the laser frequency by an amount equal to the Bohm- 
Gross frequency 




3KT e 


m 

e 



(2-31) 
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These lines are called the " satellite lines" and their contribution to 

-2 

the total scattering cross-section for a » 1 is proportional to a 
If the plasma is "cool", we see that their location tells us the plasma 
density. The actual width of the satellite lines is limited only by 
Landau damping in the Salpeter theory. However, when collisions are 

24 

permitted, the necessary corrections are treated by Taylor andComisar , 

9 25 

Dougherty and Farley , Sitenko and Gurin , and especially well by Ron, 

2 6 

Dawson and Oberman 

The right-hand term of (2-27) yields a contribution called the "ion 
feature" which is not perfectly gaussian, (having symmetric maxima away 
from zero frequency) but which has a width characteristic of the ion 
temperature. Because of the preponderant ionic mass, this width will 
usually be much smaller than the corresponding electron feature width, 
and will be harder to separate from the incident radiation at <d = 0. 

In a » 1 experiments in particular, the assumptions leading to 
(2-27) may be strongly violated if there is a very large electron-ion 
drift velocity compared to electron thermal velocity or, at microwave 

22 

illumination frequencies, if T^ » T^ . Also, as pointed out by Kegel , 
if the electron velocity distribution consists of separate cold and hot 
Maxwellians, strong enhancement of the scattering may occur, producing 
a strong central spike, and a spike in the vicinity of the plasma fre- 
quency even for a « 1. All these effects result from the behavior of 
the functions which multiply the electron and ion distributions in 

equation (2-24). S(k,a>) may be arbitrarily large at any zero of the 

23 

plasma dielectric function . 

The effects of nonzero magnetic field on S(k,cc) have also been 
7 8 11 27 28 29 

treated 9 ’ y ’ * . In general, the electron feature spectrum pre- 

dicted with no field is found to be modulated by a periodic fine structure 
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at multiples of the electron cyclotron frequency, when the sensitivity 
vector k is perpendicular to the applied magnetic field. No effect 
is expected for k parallel to the field. This modulation is a mani- 
festation of the Bernstein modes in (2-24). In our experiment, no 

effect from this source will be' seen because the spacing between the 

o 

modulation peaks is of the order 0.02 A in the azimuthal spectrum, much 
o 

less' than our 5 A resolution capability. For a, » 1, the effect of 

the magnetic field is much more complicated, depending in detail on the 

2 Q 

state of coupling between the upper hybrid and Bernstein modes . 

19 

Ichimaru discusses the effect of large-scale turbulence and inci- 
pient instabilities upon the a » 1 scattered spectrum. He has been 
able to derive an asymptotic expression, valid near a zero of the dis- 
persion relation, which predicts tremendously enhanced scattering rela- 
tive to that predicted for simple thermal fluctuations. 

References 30-33 discuss relativistic corrections to the theory 

30 

which are necessary for very hot plasmas. Brown and Rose , for example, 
show that in case Qt « 1 > a unique mapping exists between the spectral 
distribution and the electron velocity distribution. In principle, this 
relationship can always be inverted experimentally, to obtain the dis- 
tribution function f(v) from the total scattered intensity S(k). Of 
course, our plasma is nonrelativist ic . 

Computations of scattering profiles according to the Salpeter theory 

34 39 

appear, for example, in Gerry and Patrick and Williamson et al . 

In the present work, we desire only to measure the electron feature 
of the scattered spectrum, i.e. a« 1. Therefore equation (2-30) will 
adequately describe our scattering results, if we have a Maxwellian, 
nonrelativist ic plasma, since the only significant deviations from ( 2 “ 30 ) 
occur within a few plasma frequencies of the laser frequency and 

= 71 <<C 1 * (2_32) 
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3. SUMMARY OF PREVIOUS SCATTERING EXPERIMENTS 


As mentioned in Chapter 1, optical scattering from free electrons 
was first observed in the laboratory by Fiocco and Thompson in I962 1 , 
using an electron beam as the scattering medium and a normal mode ruby 
laser for illumination. In 1963* Thompson^ and Funfer et al t ^ reported 
the first successful optical scattering from a laboratory plasma almost 
simultaneously. Since then, many experiments have been reported, show- 
ing the ion feature, the electron feature and combinations of the two. 

The great majority of these experiments have involved dense plasmas 
(n^ > 10 cm giving* several orders of magnitude greater total 
scattering than is possible in our experiment. 

Among the ion feature experiments, the earliest was that of Funfer 

et al., mentioned above. This work was performed in a theta-pinch where 

1 8 ~ 3 

electron densities up to 10 cm ^ were available. However, spectral 

resolution was insufficient to resolve the central peak of the ion feature. 

The satellite lines were not seen at this time, but were observed emerging 

for 0L « 1 by the same group in ly6b^ . The experiments of DeSilva 
39 

et al. ^ in the same year was the first in which both electron and ion 

15 ~3 

features were measured in the same plasma (a hydrogen arc with n ^ 10 'em 
by changing the scattering angle from a large to a small value. Even with 
this density, stray light was about 70$ of the Thomson signal at the center 
of the electron feature. 

The satellite lines were first clearly observed by Chan and Nodwell 

4o 4i 

in 1965 , and subsequently by Ramsden and Davies . In the latter experi- 

15 -3 

ment, the plasma was a hydrogen theta-pinch with n « 2 X 10 ^ cm 

e 

and scattering angles of 13*5° and 90° gave a values of 3.0 and 0.5 
respectively. Their experiment was the first to show both features with 
good resolution. 1 
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The experiment of Rohr ; ^ first showed the detailed shape of a 

17 

satellite line, in a theta-pinch with - 5 x 10 cm ^ and = ©^ 

= 32 ev . Since a was only ~ 3*0 in this experiment , the line is not 

expected to be as narrow as in the limit a » 1, yet is was found to 

o 

be about 4 times broader 15 A) than could be explained by collisions. 
Nonuniform density in the scattering region was offered as the explana- 
tion. This experiment is also important because it is one of the first 
for which the number of electrons per Debye cube was of the order of 

3 

unity. The Salpeter theory assumes n e D e » i • 

Two other interesting ion feature experiments are those of Daehler 

44 - 4 ^ 46 

and Ribe 9 y and Kronast et al . at Garching. Daehler has shown 
scattering cross-sections several times larger than the theoretical 
value, based on known plasma parameters in the Scylla theta-pinch at 
Los Alamos. He attributes this enhancement to superthermal density 
fluctuations. Kronast et al , were probably the first to show large net 
drifts of a plasma from ion feature measurements. 

Several electron feature experiments have been done, and these will 
be mentioned only as they are relevant to our own experiment or for 
their historical interest. For an exhaustive listing, the reader is 

47 48 

referred to the papers of Kunze and of Evans and Katzenstein 

Electron feature measurements are more sensitive to stray laser 
light than are those of the ion feature, because the same order of 
scattered energy is spread over approximately \J m /m times greater 
doppler width for equal ion and electron temperatures. Accordingly, 
the ratio of stray to scattered light intensities per unit bandwidth 
at the spectrum center has usually been quite large in previous electron 
feature experiments. In the work of Funfer et al . mentioned at the 
beginning of this chapter, this ratio was about 100:1 even with 


20 



17 -3 

© = 100 ev and n = 10 cm . In the experiment of Davies and 

e e 

Ramsden^ it was 1000:1 when the electron density was 5 x lO 1 ^ cm"^ and 
the electron temperature 3 ev. 

50 

The experiment of Gerry and Rose^ in 19 66 was the first in which 

the stray laser light level was reduced to the same order of magnitude 

as the scattered signal at the spectrum center, with plasma densities 

-* 

and temperatures similar to ours. Their experiment was done in an 

argon plasma formed in a hollow cathode discharge at an electron density 
13 -3 

of only 5 X 10 cm and an electron temperature of the order of 5 ev, 
However, plasma light was always larger than the scattered signal in 
their experiment, reaching a maximum which prevented measurements at or 
beyond the strong Argon II line at 6886.6 A. As a result, their data 

show Maxwellian behavior only out to 2 thermal half-widths. Their experi- 

/ 

ment is interesting because of the great care taken in its design. 

Several electron feature measurements have been done in theta-pinches 

and other devices characterized by densities or temperatures much greater 

*52 

than ours. Only one of these, the work of Forrest, et al . in I 97 O 
is very interesting from our point of view. In this experiment, con- 
ducted by the Culham group on the Russian Tokamak T3~A machine, electron 
temperatures in the range 100 ev < < 1000 ev were found, at densities 

of order 10 J cm . Their laser energy was very similar to ours (about) 

6 joules, Q-switched) . In most cases, the stray light signal was 200 
times the true scattered signal at the spectrum center. Since, at the 

*The earlier experiment of Schwarz5 1 showed a stray laser light contribu- 
te -3 

tion only one eighth of the scattered signal for =10 cm and 
electron temperatures which must have been of order 5 ev . However, only 
electron density was measured, so that the "scattered signal" consisted 
of the entire integrated scattered spectrum. 
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spectrum center, their detector has a bandwidth designed to collect 
about the same proportion of the total scattered light as does our 
detector, their experiment may be directly compared to ours, with regard 
to stray laser light level. The result of such a comparison shows our 
stray light reduction methods to be superior to theirs by a factor of 
about 3>000. Also, measurements were made only on the blue side of the 
spectrum because of limited photomultiplier quantum efficiency in the 
red . 

In all other respects, the Tokamak scattering experiment stands 
as a model for such work. Exhaustive measurements were made of electron 
temperature versus time, versus radial position and for various plasma 
parameters. Measurements show Maxwellian behavior out to 3-5 thermal half- 
widths. Standard deviation of the data points relative to a gaussian 

is very small, of the order of 5 Of course, the Tokamak measurements 

were for only one direction of the sensitivity vector k. 

A very few experiments have also been done on shock-wave plasmas. 

r 3 s4 

Two of interest are those of Patrick^ -3 in 1965 and Martone and Segre^ 
in 1969 . In Patrick 1 s experiment, electron density behind the shock 
in a magnetic annular shock tube between 10 ^ and 10 cm , and elec- 
tron temperature in excess of 20 ev were found. His experiment showed 

4 

stray laser light levels approximately 3 X 10 greater than ours, and 
measurements were made at only two wavelengths. Martone and Segre per- 
formed their experiment in a hydrogen theta-pinch with preionization, 
and found electron temperatures of about 80 ev. 

Experiments showing the expected cyclotron harmonic modulation of 

55 

the scattered spectrum in a magnetic field have been done by Kellerer , 
Evans and Carolan^, and Kronast and Benesch^. The latter experiment 
shows clear evidence of the Bernstein modes in an OC » 1 theta-pinch 
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plasma. Kellerer and Evans et^ al . have seen the same effect for moderate 
and low a respectively. 

erg cq gQ 

The experiments of Ringler and Nodwell^ should be mentioned 

as an example of low measurements which show anomalous results. Although 
their a value is 0.4 (clearly not « l), one still expects a basically 
gaussian scattered spectrum from their steady state electric arc plasma. 
Instead, several resonance peaks were found at multiples of the plasma 
frequency and these could not be explained as due to density nonuniformities 

61 

in the scattering region. Kegel has offered the explanation that at 
least one of the observed peaks could be due to the presence of hot and 
cold Maxwellians in the electron velocity distribution. 

A very few continuous laser and plasma experiments have been done. 

62 

A first experiment of this kind was done by Farrow and Buchsbaum in 1965 

The experiment of Koons and Fiocco^*^, in 1968, is of far greater 

interest to us, because it is the first in which electron temperatures 

parallel and perpendicular to the field direction in a magnetized plasma 

were measured. Their experiment was done in a reflex discharge in argon, 

o 

with an illumination wavelength of 4880 A. Their results show a longi- 
tudinal temperature about 2.5 times larger than the transverse tempera- 
ture. Measurements were made for doppler shifts out to a maximum of 
2 thermal half-widths. Except for the type of discharge, our experiments 
are sufficiently different that a comparison of results is pointless. 

In the first place, the length of their plasma column is only twice its 
diameter. In our experiment , we have been careful to make measurements 
at a point about 15 plasma diameters distant from the reflex cathodes, 
in order to avoid the turbulent cathode-fall region near the electrodes. 
Also, Koons and Fiocco have a plasma which is only 1$ ionized so that 
it is completely dominated by electron-neutral collisions. In contrast. 


23 



the density of neutrals in our plasma is never much greater than the 
electron density. As explained in Chapter 4 , completely different 
transport mechanisms and instabilities exist for reflex discharges in 


the collisionless and collision-dominated regimes. Finally, their 
magnetic field is 1 kilogauss, whereas most of our measurements are made 
at quite low fields. Equipment scattered light is also quite large in 
their experiment, being some 20 times the scattered signal. This was 
probably due to the use of flat black rather than shiny black stray 
light trapping surfaces. 

Toward the completion of our work, we became aware of independent 

work in a continuous hollow-cathode discharge being completed at M.I.T. 

6*5 

by McCormick . This experiment is also set up to measure longitudinal 
and azimuthal temperatures in a plasma which has similar density and tem- 
perature to ours. His laser stray light level is 5 times the Thomson 
signal at spectrum center, or about 100 times the level we have achieved. 
Also, the light collection optics in this experiment are ±^ 3 - 3 ? giving 
3 times our scattering angle uncertainty and thus significantly greater 
’’cross-talk' 1 between the parallel and perpendicular measurements that 
we have obtained. To date, the M.I.T. experiment has found equal- 
temperature Maxwellians for the longitudinal and azimuthal distributions 
as we have, but measurements were made only out to about one rather than 
6 thermal half-widths. 
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4. DISCHARGE CONSTRUCTION AND PLASMA PARAMETERS 


4.1 The Penning Discharge 

The plasma which we study in this experiment is transient, highly 
ionized, moderately dense and characterized by warm electrons and cold 
ions. It is created in a pulsed Penning or Phillips Ionization Gage 
(P.I.G.) type linear discharge operating with two hot cathodes in a low 
to moderate magnetic field. 

The Penning discharge is essentially an electron beam-generated 
plasma immersed in a magnetic field parallel to the beam, but is distin- 
guished from other such discharges by the magnetron-like guns which pro- 
duce the beam. Figure 3 , which shows the construction of one of our 
cathode-anode assemblies, may be used to illustrate the basic principles 
of a Penning electron gun, although we should note that many quite 
different designs are used. 



Figure 3 

Construction of PIG Gun Assembly 


electron 

flow 


A large electric field with both radial and axial components is 
impressed across an annular anode-cathode assembly. The cathode may be 
thermionic or cold, in which case the intensity of the plasma ion 
bombardment determines the magnitude of the beam current densities. If 
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the gun is now immersed in an axial magnetic field B , electrons emitted 

z 

by the cathode in vacuum will be subject to an E X B drift with velocity 

v = cE /B (so long as B is large enough to keep the gyroradius less 
0 r z z 

than the cathode-anode spacing) and an axial acceleration a ^ = -eE^/m^. 

The result is an energetic hollow beam in the z-direction. After leaving 
the gun the vacuum trajectory of the beam electrons will follow the 
magnetic field lines and may rotate slightly due to the radial self field 
of the beam. An electrostatic reflector at the opposite end of a discharge 
tube from such a gun can cause the beam electrons to spiral back and 
forth along the tube several times before colliding with the anode. 

If now neutral gas is introduced to the discharge tube, the primary 
electrons return to the anode earlier, as the result of electron-neutral 
collisions. The advantage of this type of discharge is that the pri- 
mary mean free path is initially determined by ionizing collisions down 

-12 *66 

to neutral pressures as low as 10 torr. For this reason, the 

Penning discharge generates highly ionized plasma over a wide range of 

density. The primary limit is an upper one ^hich is due to the combined 

effects of recombination, diffusion and reduced primary range. In our 

discharge, the upper limit for good ionization efficiency in hydrogen is 

about 2| j. neutral pressure, while the discharge remains uniform to much 

higher pressures. Of course, the primary range can be many centimeters 

1 3 

in a plasma (e.g. n = 10 , T =1 ev, X ^ = 1.4 mm) whose thermal 

v e 9 e 7 mfp / 

particles have very small range, because (a) the primaries can easily 
have energies of the order 300 volts, and (b) the range in plasma is 
proportional to energy squared. 

* 

In fact, since collision frequency is proportional to density, the 
device is particularly useful for measuring low pressures. The Phillips 
Company first marketed such a gage, hence the name Phillips Ioniza- 
tion Gage or PIG, after F .M . Penning^ demonstrated the principle. 
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Hooper gives an excellent review of Penning type discharges. He 
warns that no adequate theory exists for the PIG discharge, so that the 
studies which have been done are necessarily semi-empirical. This state- 
ment can be made even more strongly for highly ionized thermionic- 
cathode PIG discharges. The reason for this is the very complicated 
state of the PIG plasma during the discharge and the associated diffi- 
culty of studying the effect of changes in isolated parameters. Strong 
instabilities and highly non-Maxwellian particle distributions may exist, 
and these can cause relatively strong electric fields even in fully 
ionized plasmas. Such electric fields, combined with ordinary density 
or temperature gradients lead to completely nonclassical particle trans- 
port modes across magnetic field lines so that the usual ideas of mobi- 
lity and diffusion velocity are useless. 


4.2 Classical and Anomalous Particle Transport 

The behavior of plasma in a magnetic field depends in an essential 
way on whether or not collisions are dominant. For purposes of dis- 
cussion, the steady-state fluid equations for electrons and ions may be 
69,70 

written 


nm|-|^- + v-7^ v = nq^E + — j - V • (nmuu) - nravv (4-1) 


where n, m, q are species density, mass and charge, v is macroscopic 
drift velocity and uu the stress dyadic, taking account of the neutral 
component in the case of partially ionized plasma or of Coulomb colli- 
sions in fully ionized plasma through the collision frequency v. If we 
now look just at situations in which the plasma is isotropic (ignoring 
viscosity effects) 7 • (nmuu) = 7p and if we ignore quadratic velocity 
terms, the simpler equation 
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dv j X B 

p at = “3— 


Vp - nmvv 


( 4 - 2 ) 


results for which a solution may be written for each of the two fluids 
in terms of their density n and the electric potential 0, in the 
form 


-v = =- ( D** V) n + (~-V)0 . 


( 4 - 3 ) 


With T = 1/v as a collision time (coulomb or neutral collisions), 
we have: 
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( 4 - 4 ) 


where D and p are the diffusion and mobility tensors, respectively, 
in a coordinate system in which magnetic field B lies along the 
third or z-axis and is assumed homogeneous. Further gd is the species 

__ Q D 

gyrof requency given by a) = . We have written this formalism in 

order to show the transition from collision-dominated to collisionless 
diffusion and mobility. This transition is accomplished by the substi- 
tution 


1 --|- ( 4 - 5 ) 

nq t 

for the resistivity of a fully ionized plasma. The diffusion coeffi- 
cients for the two cases are as follows: 
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(4-6) 


The mobilities are obtained by making the formal substitution kT -» q, 

Considerations such as these lead to the statement that classical 

diffusion across magnetic field lines in a fully ionized plasma is pro- 
2 

portional to 1/B . However, in this form it is also easy to see that 
enhanced diffusion at the rate 1/B can also occur if there are signi- 
ficant electric fields present in the plasma in the direction perpendi- 
cular to B. Relatively large electric fields (several volts/cm) can 
exist in ordinary laboratory plasmas, due to ordinary particle transport 

and instabilities, ignoring sheaths for the moment. In fact, recent 
72 

results show that such fields may be the rule rather than the excep- 
tion . 


In early plasma experiments, an attempt was made to distinguish 

between "classical" diffusion and the (1/B) -proportional diffusion which 

was first seen empirically and labeled "Bohm diffusion" . In fact, Bohm 
71 

suggested that the appropriate perpendicular diffusion coefficient 


should be D 


kTc 
l6qB ' 


4.2.1 Plasma Instabilities in the Penning Discharge 

The formal similarity between the Bohm diffusion coefficient and 

the Hall coefficient D±; — in the collisionless case shows the practical 

EXB 

importance of self-induced electric fields in a plasma. In fact, in a 


Diffusion along the magnetic field is normally not defined, but the 
corresponding mobility is. 
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cylindrical plasma, an azimuthal electric field due to charge separation 

would produce such a diffusion in the radial direction. 

7^5 74 

Such considerations led Kadomtsev and Nedospasov , Hoh , Guest 

75 

and Simon and others to develop a theory for "anomalous" diffusion 
featuring an instability known as the flute, screw, E x B or neutral 
drag instability. The model chosen was the Penning discharge in most 
cases, because of (a) the possibility of balanced longitudinal electron 
and ion currents (no net axial current), (b) the presence of a strong 
homogeneous longitudinal B field, (c) the presence of strong radial 
electric fields, and (d) the existence of much experimental data showing 
both enhanced diffusion and evidence of rotational instabilities of the 
type discussed in such devices. In such a discharge, the radial electric 
field drives both electrons and ions azimuthally at the velocity 

v = E c/B . In these early theories, a weakly ionized plasma was 

considered so that the neutral drag, that is the differing mobilities 
for electrons and ions colliding with neutrals, causes azimuthal charge 
separation when the plasma is attempting to move at this drift velocity. 
This charge separation generates an azimuthal electric field, which 
combined with a radial density gradient can grow arbitrarily large above 
a critical value of magnetic field. This electric field, of course, is 
all that is required t ( o drive large quantities of particles radially 
outward, yielding enhanced diffusion. 

If, as in our case, the degree of ionization is high, Hoh and 
Lehnert have shown that the screw instability still occurs quite readily, 
with a new charge separation mechanism due to the combined effects of 
net longitudinal current flow and collision-dominated ions. Ion colli- 
sions with the plasma tube wall resulting from large ion gyroradius can 
produce this effect. 

Other instabilities caused by finite resistivity or plasma density 
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gradients are undoubtedly present in the highly ionized Penning discharge, 
and can be responsible for the enhanced transverse particle loss rates 
that have been observed. Experimentally, the discharge at moderate and 
high densities operates in two basic modes. Below a critical magnetic 
field of the order of 200 gauss, classical or at most Bohm diffusion is 
observed. Above this field, the particle loss rate typically increases 

77 7ft 7Q 

with increasing magnetic field ’ ’ . This is not even diffusion any 

longer, but the expulsion of particles due to strong instabilities 
whose quantitative prediction has been impossible to date. 

In the region below the critical field, in which most of our 
scattering measurements are reported, the plasma is reasonably quiescent 
but may well have nonthermal distributions of electron velocities. We 
have felt that an experiment would be of interest which could determine 
whether or not the electron velocity distributions in the quiescent 
PIG discharge are Maxwellian, and whether significant anisotropy exists 
between the azimuthal and longitudinal distributions. 
k . 3 properties of our Penning Discharge 

We operate a discharge (see Figure 4) with two opposed hot cathode 
Penning electron guns at either end of a pyrex discharge tube approxi- 
mately 1 meter in length and 7*6 cm in diameter. The discharge is run 

/ 

in a pulsed fashion in hydrogen, at pressures between 0.5 /and 2 microns. 

/ 

Hydrogen was selected because of the absence of line radiation in its 
spectrum over the possible range of scattered radiation wavelengths, and 
because of the simplicity of the atom ( see Chapter 5 ) . 

In the central region of the discharge tube, a blackened stainless 
steel scattering chamber is inserted to allow connection to stray-light 
trapping chambers which will be discussed later. All metal parts are 
allowed to float electrically to avoid interference with the Langmuir 
probe. 
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magnet 



Construction of the PIG Discharge 






The electron guns are as in Figure 3> and consist of a 2 cm diameter 

resist ively heated cylindrical oxide-nickel mesh cathode of total surface 

? 

area 30 cm , surrounded by a 3 cm diameter water-cooled anode sleeve 
and backed by a tantalum disc reflector at cathode potential. Each gun 
is capable of pulse currents up to 300 amps, temperature limited. 

Gas enters at one end through a calibrated leak and is pumped 
out the other, however at such a low flow rate, about 2 liter/sec, 
that pressure differentials along the discharge tube are not signi- 
ficant . Impurity levels from virtual and real leaks remain less than 
0.01 micron-liter/sec. 

13 -3 

Typical operating parameters are: density n = 2 X 10 cm 

© 

T = 5 ev; T. = 0.5 ev; Debye length D = 5n» magnetic field strength 

B = 190 gauss; relative ionization p = 75$ » CJDT = 175 > electron mean 

' © 

free path X e = 2 cm, discharge duration = 250 (isec. With even a 
moderate magnetic field (50 gauss), the plasma forms a fairly well de- 
fined column with roughly rectangular radial density profile and outer 
radius slightly greater than the anode ring radius. 

The discharge is operated from a vacuum-tube pulser whose own 
temperature-limited output provides stabilization for the discharge 
current and a nearly rectangular current wave form into the electron 
guns. The rapid shutdown of the excitation allows afterglow measure- 
ments to begin immediately. 

Since the longitudinal resistivity of our discharge plasma is 
comparable to that of carbon, the main part of the plasma is a highly 
ionized positive column whose potential is within a few tens of volts 
of the anode potential . Most of the discharge voltage appears across 
a cathode fall region which is the main mechanism for primary electron 
injection in this operating mode. To facilitate probe measurements, 
the discharge is operated with anode grounded. 
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A typical plot of electron density and temperature vs time at 
the center of the plasma column is shown in Figure 5 . In this figure, 
time begins at the end of the discharge pulse, which has been going 
on for 200 \isec prior. These measurements were made by Langmuir probe, 
and the density measurements have been verified by microwave inter- 
ferometry. Conditions were as near as possible to those used in the 
actual Thomson scattering measurements to be reported. 


3^ 
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Figure 5 

Electron Density and Temperature vs Time at End oi Discharge Pulse 

(Langmuir Probe) 
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5. THE EXPERIMENT 


5 .1 Orientation 

Figure 6 is a plan view of the main part of the scattering appara- 
tus ; showing just the mechanical arrangement of the Penning discharge, 
magnets, light ports and scattering chamber. Figure 7 is a comple- 
mentary diagram, showing only the optical arrangement used. Both should 
be viewed in conjunction with the wide angle photograph of the completed 
experiment, Figure 8. 

The Penning discharge assembly is the structure which threads the 
magnet coils from left to right in Figure 6. Its overall length is 1.8 
meters; the plasma column itself is 1.0 meters long. 

Attached to the Thomson scattering chamber in the central region 

of the discharge are three stray light trapping baffle tubes (see section 

5 . 3 ), part of the vacuum system, at angles of 60 ° to the plasma axis. 

These convey the laser beam through the scattering region in either of 

two directions k._ or k. . and the scattered light to its detection 
ll r2 

system. The short section on the fourth flange of the scattering chamber 
is a viewing dump, which is a "black hole" designed to provide a proper 
background against which to view scattered light . The beam focussing 
and viewing lenses (L^> and in Figure 7) are located just outside 
the windows of the three baffle tubes. The top of our Langmuir probe 
is seen in the center of the scattering chamber. 

The scattering angle between incident and scattered wave directions 

s o o 

is seen to be either 60 or 120 , depending upon which illumination direc- 
tion is used* This choice of angles is primarily dictated by mechanical 
constraints (permissible field coil locations). Actually, any combina- 
tion of angles such that 8^ + with the discharge axis bisecting 

the angle between illumination and scattering axes would serve to give the 
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Figure 6 

Schematic View of Scattering Experiment from Overhead (Optics Removed) 
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Schematic View of Illumination and Detection Optics 



Figure 8 

Wide Angle Photograph of Experiment 
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necessary sensitivity directions and for measuring azimuthal 

and longitudinal velocity distributions (see Figure 2). The parameter 
a in the two directions is, respectively, of the order 0.03 to 0.02 
during the discharge. 

We could equally well have used a fixed illumination direction and 
two observation ports. The advantage of doing so is measurement 
simultaneity. The disadvantages are the cost of duplicating the scattering 
detection system, including spectrometer and photomultiplier, and, more 
important, the fact that relative errors due to miscalibration of the 
detection systems would add rather than cancel. In this experiment, both 
the laser and tne plasma are highly reproducible from shot to snot over 
long periods of operation, so we may use the single port method. 


5 .2 The Laser Light Source 

The light source for our scattering experiment must meet several 
requirements, some of which are unusual for commercial equipment . The 
most obvious of these are sufficient energy and spectral energy bright- 
ness (energy per unit wavelength, area and solid angle). By rewriting 
equation (2-5) in terms of the differential Thomson cross-section and 
scattering length £, we see that the minimum illumination energy re- 
quired to produce at least N scattered photons in detector solid 

s c 

angle dQ, and bandwidth do) at doppler frequency co is 


L (co) = 
o' y 


hv N 
o sc 



joules . 


( 5 - 1 ) 


More specifically, if we assume a Maxwellian, nonrelat ivistic plasma, 
we can use equations (2-20) and (2-27) to calculate the cross-section 
and the required illumination energy at the center of the electron feature 
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spectrum in the scattering plane (a « 1, CD = 0, sinY = l). This gives: 


L 

oe 



CD 

— joules . 


(5-2) 


Since, in any experiment of this type, *the detector bandwidth will be 

determined with resolution of the thermal doppler width of the scattered 

feature in mind, we replace (cD^/dcn) by a constant R indicating the 

resolving power. Then, in practical units for this problem (wavelength 

12 _3 

in microns, density in 10 cm , energy in joules) we find: 


L = 4-. 4 X 10 
oe 
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R N 


sc 


^Vel2^ 


joules 


(5-2a) 


for the minimum electron feature energy, 


CD. 


The same procedure for the ion feature (oc » l) with R = - gives: 

x dcD 


L . =E t j 1 + z — ^ j L 
oi Z V T^ J oe 


(5-2b) 


Stated in this way, the single number N gc implicitly includes any nec- 
essary calculations of transmission loss and signal-to-noise ratio in 
the detection system. 

As explained in Chapter 6, the only noise we have to deal with is 

statistical noise due to the discreteness of the tiny electrical signal 

present in the first stage of light detection. However, even with 

multiple measurements under the most favorable noise conditions, a 

practical minimum for this signal at the 1/e point of the electron 

feature spectrum is 10 electrons, or 27 at the spectrum center. Assuming 

a probability of 0.07 tor conversion of a photon into a photoelectron 

at the detector and an average optical transmission of 0.20 for the 

detector optics within the pass band we require N = 1,930 photons. 
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The typical parameters for the experiment we plan are n = 20, 

= 0.694, dO = 0.01 sr, R = 10 and <£, = 0.8 cm, so that the required 
illumination energy is 0.77 joules. 

The foregoing is precisely the reverse of our actual design pro- 
cedure, but equation ( 5 - 2 ) provides a particularly useful way of looking 
at the energy requirements of a scattering experiment. 

While an ordinary light source could supply this amount of energy, 
only a laser source has the necessary brightness and monochromaticity. 

The monochromaticity requirement arises from the fact that the illumi- 
nation energy must all lie in a band smaller than the spectral resolution 
desired. Usually, the most difficult part of this limitation is 
achieving an absolute stability of the illumination spectrum center wave- 
length commensurate with the permitted bandwidth. Brightness is needed 
because the beam solid angle and focal spot size both must be small 
enough to satisfy obvious bounds on incident beam direction uncertainty, 
and spacial resolution in the scattering volume. Further, the bright- 
ness of the illuminated plasma should exceed its natural brightness 
within the scattered spectrum. As an example, the spectral energy 

g 

brightness of the laser used in our experiment is about 4 X 10 joules/ 

2 o 
sr-cm -A. 

Up to this point, we have deliberately avoided mentioning power or 
time duration. In our experiment, rapidly varying plasma parameters re- 
quire a measurement interval of 1 microsecond at most. In any case, it 
is seldom desirable to obtain the required illumination energy over a 
long measurement interval at low power because of the complicated 
apparatus required to elicit the scattering signal when it is deeply 

** 

The resolution needed is a few Angstroms in an electron feature experi- 

o 

ment, or a few tenths A in an ion feature experiment. 
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buried in noise. At present, commercial lasers designed for continuous 
operation feature power levels six orders of magnitude below that of 
Q-switched lasers, requiring plasma stability over a correspondingly 
longer period. In short, nothing can match the joy that comes from 
having a real signal-to-noise ratio of 20 db! 

The source of illumination for our scattering experiment is a Korad 

model Kl-Q Q-switched ruby laser, specially modified for spectral purity. 

The laser cavity was set up to produce a 25 to 30 nanosecond, 1.5 joule 

flash of light at 69^3 A, about 50 % of the maximum available Q-switched 

energy. Typical intensity in the O .95 cm diameter output beam is 60 

2 

megawatts/cm . With these operating parameters, the laser can be fired 
ten times a minute, and the critical cavity parts have a life expectancy 
of about 50,000 shots . 

Other characteristics of the laser output are: beam divergence, 3*0 

milliradians for 95 $ beam energy; center wavelength uncertainty and 
spectral width, 0.03 A or 2 Ghz, greater than 95 $ linear polarization 
and 3 $ standard deviation of the light pulse energy from shot to shot 

at 10 shots per minute. The resulting source spectral brightness is 

l4 2 o 

about 1.5 X 10 watts/sr-cm -A. 

The low bandwidth was achieved by the use of commercial resonant 
Fabry-Perot etalons for front and rear cavity reflectors. The reflecti- 
vity of each etalon is a comblike function of wavelength, having several 
maxima and minima within the wavelength interval for appreciable gain 
in the ruby rod, but the spacing between the teeth is slightly different 
for the two combs. By tuning one etalon relative to the other, their 
reflectivities can be made to reinforce significantly in a single narrow 
cavity gain maximum corresponding to only five or six longitudinal modes. 
In our laser, the rear etalon is fixed, while the output etalon may be 
temperature tuned. Correct adjustment of the etalons for narrowband 

^3 



Figure 9a 

Laser Beam Pattern, z = 3 meters 



Figure 9b 

Laser Beam Pattern, z = 5 meters 
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operation was verified with a spectrum analyzer, using measurement ranges 

o o o 

o! 1.20 A, 0.24 A, and 0.06 A. 

The laser’s excellent output energy stability is due mainly to 
accurate control of lamp input energy through the use of oil-filled 
(rather than electrolytic) energy storage capacitors and well-designed 
capacitor charging circuitry. In addition, the flashlamp and ruby rod 
are water cooled, with ± 1°F regulation. This is important for stability 
of output power, energy and wavelength and for low beam divergence. Un- 
controlled heating during the flashlamp pulse can easily cause output 

o 

wavelength shifts of 1 or 2 A in ruby and, with the natural temperature 
dependence of ruby gain plus the wavelength-dependent cavity gain with 
mode selection, very erratic output power. Heat also causes significant 
changes in beam divergence due to warping of the ruby. The same cooling 
system is used to control the output etalon temperature. 

Another advantage of the laser head design for our purposes is 
the use of a helical flashlamp and diffuse lamp cavity reflectors, which 
results in unusually uniform pumping of th’e ruby. This is demonstrated 
in Figure 9 which shows actual photographs of the laser beam pattern on 
ground glass, after it has traveled 3 and 5 meters, respectively. Even 
at these distances, "the beam diameter shows no appreciable increase. The 
beam expander used to achieve this degree of collimation is discussed in 
the following section. \ ; - 

5-3 Plasma Illumination 'Optics 

The laser light output is sent through a 5^ beam-expanding collimating 
telescope, in Figure 7 - The function of this device is to reduce 

both the beam divergence and power -density of the beam, while eliminating 
laser pump light and other stray radiation outside a narrow forward cone. 
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The telescope has an effective f-number of 10.0, chosen to permit 

rapid collimation with simple lenses. Lens focuses the laser energy 

on a 0.55 mm diameter spot in the center of D^, which is a sharp-edged 

O .95 mm diameter pinhole in 0.005" tantalum sheet. The power density 

2 

at the focus is of the order of 40 gigawatts/cm , so the pinhole was 
deliberately made too large to absorb a significant amount of the main 
beam energy. Instead, all radiation in a cone of half-angle 4.9 milli- 
radians (roughly 1/4 ) about the forward direction is transmitted. The 
laser beam waist in the plasma is then determined by considerations other 
than the size of 
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As pointed out by Evans and Katzenstein , such focal power densities 
are sufficient to cause a spark in air, which would absorb and scatter 
a large and unpredictable amount of the beam energy. This is cited in 
their paper as being a basic difficulty with fast collimation, but we 
have found that evacuating the focal region to 'about 5 microns solves the 
problem. 

After expansion and collimation, the beam has a diameter of 4.75 cm > 

the divergence has been reduced to 0.6 milliradians for 95 $ of beam energy, 

2 

and the average power density is only 3 megawatts/cm . Overall trans- 
mission of the collimator for the raw laser beam is 69 $. This power 
level has proved quite safe for the multilayer dielectric coatings used 
for mirror and antireflection surfaces in the laser beam optics. Further, 
the very low expanded beam divergence allows us to work with an essen- 
tially constant diameter beam throughout the experiment, which greatly 
simplifies handling of stray laser light. 

The collimated laser beam next passes to the double-sided mirror 
M which can be rotated to direct the beam through the plasma scattering 
region either along the path shown in Figure 7? or in the reverse direction. 


46 



Because of the symmetric train of optics L the beam is 

focussed on the same volume of plasma for either incidence direction, 
to within a maximum misregistration of ± 0.5 mm . In the scattering 
region, the laser beam is vertically polarized. The beam waist diameter 
in the focal region is intended to be 1.0 mm for 95 $ of the beam energy, 
and has been deduced f rom measurement s of the Rayleigh scattering inten- 
sity profile at the focus to be 0.95 ^ 0.05 mm. The resulting energy 

2 

density in the scattering region is then, typically, I 5 O joules/cm . 

The amount of plasma heating expected at this energy density will be 
discussed in Chapter 6. 

It should be mentioned here that the effective f-number of the 
focussing lenses is completely determined by the desired beam waist 
diameter for a beam of known etendue (the product of illuminated solid 
angle and beam cross-sectional area). For example, our laser has an 
inherent etendue = 1 . 9 ^- X 10 y sr-cm which gives an empirical 

_ o 

relationship d = 5.62 X 10 J f cm for the waist diameter. This was 

03 f 

the primary reason for picking an f-number of 15 for these lenses, although 
this choice incidentally results in a quite acceptable contribution of 
2° to the standard deviation of the scattering angle. The one degree of 
freedom which exists is the expansion ratio of the beam collimator, 
which, when chosen determines the focal length of and L^. Our ratio 

of 5-0 results in a ’Jl cm focal length, which is a desirable figure from 
scattered light considerations. 

In order to scan the scattering region across the plasma column, 
we chose a very direct method in which mutually calibrated, parallel-axis 
micrometer drives allow simultaneous vertical translation of the mirror- 
lens assemblies and M^,L^ for the laser beam, and Mg,L,_ for 

the detected beam. A total translation of 3*8 cm is permitted. 





All lenses and windows in the laser optical train are made from 
objective quality "Vycor 11 glass for thermal stability and uniform refrac- 
tive index and have a surface figure accurate to within one wavelength 
of sodium M D M light. They also have a resonant multilayer dielectric 
coating which reduces the maximum reflectance at each surface to 0.1$ 
at the laser wavelength. All mirrors are borosilicate glass flat to 
1 wave per inch and provided with a resonant multilayer dielectric re- 
flecting coating. For light at the laser wavelength, incident at 45° to the 
surface normal, these mirrors provide a reflectivity of at least 99.9 $ (see 
Figure 10). However, they are relatively transparent to light of substan- 
tially different wavelength, and so enhance rejection of laser pump light.. 

As a result, at least 98$ of the laser light output at the colli- 
mator succeeds in passing through the experiment and bouncing off the 
second surface of into the beam dump. 

The dump is a wedge constructed of two plates of Schott BGl8 red 

absorbing glass, one piece intersecting the beam at the Brewster angle 

-5 

and the other parallel to and outside the beam. About 5 X 10 of the 
incident light is transmitted to the beam monitor photodiode by the 
Brewster plate, giving an electrical signal which is integrated to deter- 
mine the relative energy of each laser pulse. 

This is a somewhat unusual approach in that the laser light is 
absorbed far from the scattering region, and in such a way that return 
scattered light from the dump is effectively eliminated by geometry. 

The remaining sources of equipment scattered light at the laser wavelength 


are then the mirror-lens pairs 


m 3 l 3 


and the adjacent scattering chamber 


windows, shown in Figure 11. 
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Figure 10 

Typical Mirror Transmissivity vs Incidence Angle at 69^3 
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5.4 Stray Laser Light 


The scattered light spectrum in this experiment has a total inten- 

-14 

sity of the order of 10 of the illumination intensity because of the 
modest electron density and the small range of scattering angles we are 
permitted to use. For this reason a source of noise at the laser wave- 
length is laser light diffuse-scattered into 2jt sterradians by the 
optical elements closest to the plasma and again scattered into the detec 

tion optics by the equipment itself. Even if the spectrometer has a 

* 5 

rejection ratio of 10 , it cannot adequately discriminate against such 
light that has not already been reduced to 10~^ or so of the incident 
light level.' 

Figure 11 shows the internal construction of the stray light baffle 
tubes which accomplish a major part of the trapping. A similar but 
shorter baffle tube is used at the scattered beam port (see Figure 6). 

In essence, these tubes are evacuated, rough-collimating light ports 
with perfectly black walls.. The trick is to make the walls "black", at 

least for red light. This is accomplished by a telescoping series of 

/ 

thin, shiny black irises and light-tight spacers within each tube. The 
aperture of each iris is oblong and just slightly larger than the local 
cross-section of the collimated beam, allowing for beam translation. The 
frequency of irises along the beam axis is determined geometrically by 
the requirement that any uncollimated ray originating within the beam 
volume at the entrance window, or at points more distant from the beam 
focus, will be trapped by multiple reflections between adjacent irises. 

* 

The "rejection ratio" of a band-pass filter is defined as the ratio of 
light transmitted at the point of peak transmission to that transmitted 
at adjacent frequencies where the filter is supposed to have zero 
transmission. For grating monochroma ters, this ratio is seldom greater 
than 1000. 
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Figure 11 

Internal Construction of Stray Light Trapping Tubes 


Typical rays will be attenuated by 6 or 8 orders of magnitude before 
escaping from the wall region. 

This technique depends upon the use of specular rather than diffuse 
iris and spacer surfaces. An ordinary "flat black" surface would scatter 
a significant portion of the stray light back into the focal region 
at the first encounter. Of course, the entrance windows must be tilted 
slightly out of the plane perpendicular to the tube axis, so that the 
small amount of properly collimated light specularly reflected from 
their surfaces may also be trapped. 

Because the baffle tubes communicate directly with the Penning dis- 
charge, considerable effort was necessary to find a shiny black coating 
for their internal parts which would not outgas. We have developed a 
sequence of surface treatment procedures which permit the use of commer- 
cial paints for this purpose, achieving an outgassing rate of 0.005 

micron- 1 iter/ sec per square meter of surface at room temperature. System 
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partial pressure due to paint is about , 3 X 10 torr. In addition, a 

very small amount of auxiliary pumping is done on the baffle tubes to 

maintain a net gas flow away from the scattering region. 

Beyond merely trapping stray light, the entire experiment has been 

designed to reduce stray light flux by taking advantage of the inverse 

square law wherever possible. We have mentioned the use of a distant 

laser beam dump. Similarly, light from the laser head must travel appro- 

o 

ximately 3*5 meters through 5 right-angle turns and a 12.5 deviating 
prism (not shown in Figure J) before entering the laser beam port. 

Probably the most important measure of this type has proved to be 
locating the laser beam port windows 70 cm from the scattering region. 

The aperture next to the plasma in each of these ports then intercepts 

-4 

only 1.3 X 10 of the light scattered from a point in the entrance pane. 
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It might appear that nothing could be gained in this way, since we are 
dealing with fixed beam f-numbers for which the illuminated window area 
increases as the square of the distance from the focus. However, it 
should be remembered that the scattered intensity at the window is pro- 
portional to the incident intensity, and that both decrease as the dis- 
tance squared, leaving a total scattering from the window which is indepen 
dent of its location. 

80 

In fact, if we are to believe the military specifications for 
optics, which treat only gross imperfections such as scratches and bubbles 
a piece of given quality should have a per unit area scattering cross- 
section which actually varies inversely with the diameter of the piece, 
leading to an even greater stray light reduction! 

The final surface which must receive careful attention to achieve 
a low stray light level is the background against which the scattering 
is viewed. If this surface (the "viewing dump" ) returns no light at all 
to the detector, we have achieved our purpose. The original model for 
a viewing dump is called a Rayleigh horn, after its inventor, and is just 
a hollow cone of black glass, curved into a horn-shape so that the point 
cannot be seen when looking at the inside. Its blackness depends upon 
the absorption of stray light in multiple specular reflections between 
highly polished black surfaces in much the same way as did light trapping 
in the baffle tubes. 

However, a Rayleigh horn would be large and awkward in this experi- 
ment, so a simple wedge like the laser beam dump is used, but with a 
smaller dihedral angle of 20° between the Schott BG18 glass plates. As 
illustrated in Figure 12, 8 or 9 bounces are required for stray light to 
enter the viewing beam. The line of intersection of the wedge plates 
is out of sight of the detector. Further, the top and bottom of this 
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wedge are covered by glass plates, resulting in a sort of flattened 
Rayleigh horn which is tailored to the oblong footprint of the scanned 
viewing beam. Finally, the dump is backed by a black surface of the 
same shape, and the whole structure placed at the rear of a short, dead- 
end port provided for the purpose. 

The last element of our stray-light reduction program is the elimi- 
nation of light transmitted to the scattering chamber by total internal 
reflection in the glass sections of the Penning discharge tube. This 
was done by painting the outside of these sections with a commercial 

black paint which reduces the reflectance at the exterior surface to 

-6 

about 10 by accurately matching the glass refractive index. 

The ultimate test of these techniques is the measured equipment - 
scattered light level at the laser wavelength. This light level is one 
of the lowest attained in any experiment of this kind, being equivalent 
to Rayleigh scattering from no more than 3 microns of neutral nitrogen 
gas, or about 20^ of the typical center-spectrum Thomson signal. The 
actual value in a particular experiment can be 2 or 3 times this figure 
if the entrance beam optics have been allowed to become dirty. The main 
point we wish to make is that ours is the only low-density experiment 
where accurate Thomson me asurements have been possible at the laser wave - 
length. As discussed in Chapter 3j our equipment scattered light level 
is about 3 jQQQ times better than that reported by the Culham group in 
their Tokamak measurements. 


5-5 Scattered Light Detection Optics 

The Thomson scattering volume is- defined by the intersection of 
the volume illuminated by the focussed laser beam and the volume observable 
from a point on the detector surface. As shown in Figure 7* the observable 
volume is determined by the image of aperture projected* on the 
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scattering region by the lens pair L^,L^. This aperture is rectangu- 
lar, with fixed width and adjustable length, and oriented so that only 
the illuminated plasma is seen, while the length of the plasma sample 
can be varied to optimize received signal strength and resolution. Lens 
L^ then returns the scattered beam to a parallel state as required by 

„ J 

the multilayer filter F^, which is the detected light spectrometer. 

Light transmitted by F^ is detected by a photomultiplier. The charac- 
teristics of this filter are discussed separately in section 5 .6. 

Aperture is a large, adjustable, circular iris which determines 

the f-number of the viewing system. It is usually opened to the maximum 
useful aperture of the system, which is an effective f-number of 8.33* 
Since the lens L^ has a 47 cm focal length, the scattered light viewing 


solid angle is at most 


dfi = 0.011 sterradians 
s 


The maximum opening of consistent with useful spacial resolution in 

the plasma and beam parallelism at is k mm x O .76 mm, giving an 

image in the plasma which corresponds to a cylindrical Thomson scattering 

volume of projected length 7 .^ mm, or true length 8.6 mm, and diameter 

8l 

0.95 mm. Therefore, the maximum available viewing etendue u in this ex- 
periment (the product of beam area and solid angle) is 


A -ll O 

$ s dfl i d = 7.8 X 10 sr-cm 


s is s 


Finally, the total standard deviation of scattering angle due to 

both the illumination and detection beam ray bundles is 3*9° or 6 . 5 $ of 

the mihimum angle used for scattering. This figure corresponds to an 

o 

uncertainty of roughly ±1 in the direction of the sensitivity vector 

_ « 

k for either of the two scattering angles. This means that there will 
be no more than a 2$ cross-talk between measurements of the transverse 


56 



and longitudinal velocity distributions. 

The rotation axis of is parallel to the long dimension of slit 

and perpendicular to the direction of linear polarization of the 

scattered light as it falls on the filter. 

82 

The Gunter bounce technique for enhancing photomultiplier sen- 
sitivity is not employed in our detection system because the effective 
gain enhancement is a quite critical function of incident light angle. 
Also, the maximum gain of 2.0 which could be had in this way was not 
worth the effort. 

Careful attention has been paid to the effects of diffraction and 
of spherical aberration in the design of all optical elements in this 
experiment . 

Diffraction plays no significant role. By design, the collimated 
laser beam illuminates no apertures. Diffracted stray laser light will 
not be able to produce measurable Thomson scattering because of its low 
intensity, and is simply removed along with the rest of the equipment 
scattered light as discussed above. Where apertures are illuminated, as 
in the detection optics, they will be lens pupils, and we might ask 
whether diffraction might interfere with image formation. However, a 
quick calculation verifies that the beam etendues used throughout the 
system are very much larger than the pupil diffraction etendues . , As 
shown, for example, in reference 83 the diffracted intensity distribu- 
tion in the focal plane of a thin lens uniformly illuminated by a mono- 
chromatic plane wave with wavelength \ and intensity I is 


liil 

1 


o 



(5-3) 


where r 




f is the lens focal length, and f^ = f/D its 


effective f-number as determined by a finite circular pupil of diameter 
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r 


D in the lens plane. This is the so-called "Airy pattern" of the lens. 

The central zone of this pattern, bounded by the first zero of the 

Bessel function, contains 84$ of the diffracted energy, and has a radius 
* 

r 1 = 1.22 f^X Q . The fraction of the total energy which is bounded by 
the nth zero of is 


I(r n ) 2 

= 1 - J 2 (r /r ) 
o v n o' 




(5-4) 


so that 95.2$ of the energy lies within the fourth zero, at radius r^ = 

4.23 fj/X . This latter radius corresponds to an etendue 
t o 


$ 


d 



(^.23)\ 2 


2.1 X 10 


-7 


2 

sr-cm 


which is independent of f^. This figure is 100 times smaller than the 
smallest etendue employed in the system. 

Sperhical aberration is not so easily dismissed. In order to keep 
the cost of our lenses below $1,500, aspheric surfaces could not be 
used. Instead, plano-convex surfaces were used throughout. Insofar as 
shape is a factor, this choice nearly minimizes spherical aberration 
when a simple lens is used to focus a parallel beam, provided that the 
beam impinges on the convex surface. For such a lens, the angular aberra 
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tion of the edge rays is 


d0 


_ n^-2n 2 +2 
ab l6n(n-l)^f^ 


(5-5) 


* 

This is the origin of the so-called "Rayleigh limit" for diffraction- 
limited optics. 
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which, for ordinary materials and red light, amounts to 


de 


ab 


145 

f# 3 


milliradians (mr) 


(5-6) 


A survey of the experiment^ in Table 5-1 shows that spherical aberra- 
tion is within reasonable limits, particularly as the typical ray will 
experience much less aberration than an edge ray. As a counter example, 
an f-number of 3-0 for the viewing lens would cause 1 6 milliradians of 
aberrations . 


Table 5 -1 

Survey of Spherical Aberration 

location 

f-number 

edge-ray spher- 
ical aberration 
(mr.) 

limiting 

divergence 

(mr.) 

nature of 
limit 

laser colli- 
mator output 
lens 

10.0 

0.15 

0.6 

design value of 
expanded beam 
divergence 

plasma illumi- 
nation lens 

15.0 

0.04 

0.6 

same as above 

plasma viewing 
lens 

8.33 

0.25 

2.0 

degradation of 
filter (Fj_) 
bandwidth 


The plane windows through which the converging illumination and 
scattered beams must pass also cause lateral aberration, but in a direc- 
tion opposite to spherical aberration, by increasing rather than de- 
creasing the effective focal distance of the edge rays. This effect has 
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also been calculated. Its magnitude is such as to eliminate aberra- 
tion in the illumination lens, and to give net viewing lens aberra- 
tion of the same magnitude listed in Table 5~1 but opposite sign, for 
the edge rays . 

5-5 .1 The Multilayer Dielectric Filter-Instrument al, Properties 

We use multilayer dielectric filters to scan the scattered spec- 
trum because they are cheap, and because they can provide 2 or 3 orders 

of magnitude greater etendue than a monochromator when a modest reso- 

o . 

lution (d\!> 1 A) and tuning range (AX/X q ~ 10$) are sufficient. 

Such a filter consists of two or more stacks of plane parallel 
quarter-wave layers of alternating high and low refractive index, with 
integral half-wave high or low index spacers between adjacent stacks. 
Their gross behavior can best be introduced by considering light trans- 
mission through a single half-wave layer with refractive index n 

K 
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In the diagram, a plane wave of wavelength X ^ \ Q is incident 
from air on the layer at angle 0 . The layer has phase thickness 

x ° 

o 

n^t =: — . The transmitted wave will be a maximum when the rightward 
moving wavefronts resulting from multiple internal reflections are 
all in phase. For example, the amplitudes at points 1 and 2 will have 
a phase distance difference 


V 


n( d-d ) = w (1 + sin a) = X cos 0 
k v 2 1 cos 0, o k 


since a = |[ 


- 20 k * If we require that n^d^-d^) = 


. cos 0. 

X k 


so that a simple tuning formula results for the first transmission 
maximum ( p = l) : 



1 


• - 4/2 

sm 0 

o 



(5-7) 


or 




using Snell's law to relate internal and external incidence angles. This 
shift toward the blue is somewhat counter-intuitive since it might at 
first seem that increasing incidence angle, and the resulting increase 
in ray path length in the film, would lead to a red shift. 

In fact, it can be shown (see, for example, reference 85 ) that the 
tuning formula for a complete filter including the quarter-wave stacks 
is also given by (5~7) to within a very small error, if n is replaced 
by an "effective filter index" n f , which depends upon details of the 
filter construction. For practical purposes, the simpler tuning formula 
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AX = -b 0 q 2 A (5-8) 

will be adequate. It is accurate to within 2$ for angles less than 20° 

o 

with our filter, using b = 1,350 A for n f = The exact formula, 

derived in a paper now in preparation, shows a very slight dependence 
on the state of incident light polarization. For our filter, the varia- 
tion amounts to 2.5$ of the tuning range at an angle of 30° • 

By differentiating equation ( 5 - 8 ), it can be seen that tight limits 
exist on the allowed nonparallelism of the incident light if the effec- 
tive filter resolution (reciprocal bandwidth) is not to be degraded. 
However, since those rays which diverge from the incident beam axis in 
the plane of incidence will generally have a much more serious effect 
on the filter resolution than those which diverge perpendicularly to 

this plane, simply writing -5\ = 2b0 60 is not sufficient. 

o o 

Two cases may be distinguished. The first, 60 q » 0 q describes 
perpendicular beam incidence. Here, 

a \ -^ (60 o )2 ° A ^-9) 

0 =0 
O 

gives the standard deviation of peak transmission wavelength for a beam 

whose rays diverge uniformaly in all directions up to half -angle 60 q 

2 

from the beam axis. If the filter has face area A cm , the et endue for 
normal incidence is then 

$ f (0) = JtA(60 Q 2 ) = a. sr-cm 2 . (5" 1 0) 

max 

In the second case, 60 « 0 , it will be useful for further dis- 

o o 

cussion to express the deviation of incidence angle 60 in terms of 
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its angular components 60^ and 80 with respect to the plane of 
incidence . This gives: 


60 ( 60 ) 2 + ( 60 ) 2 
,_1 + —Jzl _Ji_ 

0 0 2 

oo 20 


60 




In the traditional method of illuminating filters in experiments such 
as ours. 60 — 60 with maximum half-angle 60 as before, so that 

-L II O ' 


b 0 60 A 

o o 


0 ^0 
o 


and 


x / A y \ It A 2 2. 

*i (iX) - bTW a x sr ' om 


(5-12) 


This is a gross underestimate of the available filter etendue , as will 
be shown. By equating et endues , we find that the plasma sample area 

V 

we are allowed to view is only 


A = 
P 


4Af # 2g x 2 

b|'AX 


cm 


(5-13) 


if 60 = 60 . In our experiment, if we require a, = 1 A, a maximum 

± II ’ . * 

tuning range | A\| = 200 A and viewing lens f^ = 8.33; the projected 

2 

area of the plasma sample volume could be at most 2 mm . 

We are most interested in the situation where most of the divergence 

* 

occurs perpendicular to the plane of beam incidence on the filter. 

We see from (5-11) that the contributions of 60 and 60 to 60 are 
v ' 1 II 


The fT plane of incidence 11 is always defined as that plane which contains 
the incident beam symmetry axis and the surface normal. 
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o 

comparable only when 60 » 60 that is. for ( 60 ) 0 60 . 

± II ' i o || 

Our contribution has been simply to align the filter rotation axis 

in such a way that the larger of the two components of 60 will be 

perpendicular to the filter incidence plane. Then, since £» w, we 

have 60 = ^/2f and 60 = w/2f. Now, £ may be made as large as 

i II 


"t 2 = jtwf 0 

o 


(5-14) 


without further degrading the filter bandwidth. Expressed in terms of 

peak wavelength standard deviation, this degradation is now only cr ^ 

A 

O 


b0 60 A so that we have as a final result for the etendue achievable 
o h 

with the new method of illumination, 


* f '(AX) 


« A / v3/2 2 

(ct . ) sr-cm 




(5-15) 


This represents an improvement of ^IaXi/c., over the usual experimen- 

A 

O 

tal etendues, as expressed by (5-12) above. Again picking o\ = 1 A and 

A 

11° o 2 

AX = 200 A for a numerical estimate, we can now see at least 28 mm 
1 'max 

of plasma area. The actual magnitude of our filters etendue at |ax| 
o -3 2 

= 200 A is 3 x 10 sr-cm , which is a figure about 100 times that of 
a good monochromator, but only a factor of 4 greater than we require in 
our experiment. At normal incidence, the filter may easily be 1000 times 
better than the monochromator and this advantage is indpendent of reso- 
lution. Unfortunately, the filter must be tuned. 

Multilayer filters also have 2 or 3 orders of magnitude better re- 
jection ratio for light outside the pass band than do monochromators. 

This is a severe disadvantage for monochromators, since they must be 
used in series with another dispersing instrument for successful Thomson 


scattering work. 
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5.5*2 The Multilayer Dielectric Filter - Measured Properties 

The useful tuning range of a multilayer filter is limited by degra- 
dation of the transmissivity and/or resolution at large incidence angles 
Further, as mentioned in the previous section, the resolution or 
(bandwidth) 1 will be further degraded if the incident light is non- 
parallel. Of specific interest are the total or integrated transmissi- 
vity and the effective resolution as functions of illumination incidence 
angle, divergence angle and state of polarization. The integrated 
transmissivity must be known precisely so that its variation may be 
removed from the Thomson data. 

Accordingly, a recording spectrophotometer was assembled to measure 
these parameters (see Figure 14) . The parallel white light source is 
a uniformly illuminated pinhole in the focal plane of a high quality 
collimating lens. The beam divergence is controlled by the size of the 
pinhole, and may be 0 1 , -55, 1*70 or 4.85 milliradians . The smallest pin- 
hole used was 38 \i in diameter. A series of 20 spectrograms at varying 
incidence angles from 0 to 30° was repeated 6 times for each of the 4 
filters tested, in order to have measurements for the 3 beam divergence 
settings in both parallel and perpendicularly polarized light. Part 
of such a series for the filter used in our experiment is reproduced in 
Figure 15 , normalized to the incident spectral intensity Of course, 
the original spectrograms were taken in much greater detail. In this 
figure, the degradation of filter resolution at large angles is clearly 
seen, but the variation of the area under each curve (the integrated 
transmissivity) is not yet clear. 

Separate measurements of this quantity were achieved by completely 
opening the spectrometer exit slit to accomodate the entire filter 
spectrum and repeating measurements of transmitted intensity versus 
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Spectrophotometer Setup 
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Typical Filter Spectrograms at Various Incidence Angles 


incidence angle. After again normalizing to incident spectral inten- 
sity, a plot such as reproduced in Figure 16 results. 

From the filter spectrograms, additional plots were made of peak 
transmissivity and bandwidth as a function of wavelength, and of peak 
transmission wavelength shift versus incidence angle, all with state 

of polarization and beam divergence as parameters. The bandwidth and 

o 

tuning curves for the 5 A width filter used in our experiment are shown 
in Figures 17 and 18. 

Several conclusions can be drawn from all this data. 

(1) We can see that the integrated transmissivity of this filter 

is much more nearly constant than is its peak transmissivity, so that 

only a small relative correction to the Thomson data will be required 

o 

over the doppler range of about 200 A employed. 

(2) There is a significant dependence of the bandwidth and inte- 
grated transmissivity on the state of polarization of incident light. 

(3) The degradation of resolution with filter tuning is tolerable 
because it occurs in the slowly varying wings of the scattered spectrum. 

(4) The amount of degradation with filter tuning is much smaller 
if the incident light is linearly poLarized parallel or perpendicular 
to the incidence plane, rather than circularly polarized, as in most 
published measurements. This is because the splitting seen in Figure 18 
for large angles is added to the apparent bandwidth for circularly 
polarized light. 

( 5 ) The integrated transmissivity is insensitive to incident beam 

* 

divergence up to about 5 milliradians . In fact, selective tests show 
this to be true even at 11 mr. If this were not true, a different 
calibration would be necessary for every change in plasma sample volume 
dimensions. In our experiment, 60 is at most 0.8 mr . 
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Integrated Transmissivity vs Wavelength, Filter 
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Measured Tuning Curve for Filter 


Fortunately, Thomson scattered light is linearly polarized, and it 

is a simple matter to orient the filter rotation axis so that the filter 

sees only one state of polarization. (This is the same orientation that 

maximizes etendue when > w.) However, the integrated transmissivity 

of a filter for any state of polarization can be accurately synthesized 

from the parallel and perpendicular data, providing the filter behaves 

isotropically. Figure 19 demonstrates that this is true. Here we plot 

the calculated response of the filter for a polarization vector inclined 
. o 

at 60 to the plane of incidence next to the measured integrated trans- 
missivity for this situation. 

Finally, we note that the integrated transmissivity of a general 

multilayer filter is much less well behaved as the filter is tuned than 

was true for the one discussed in this section. For example, a nominal 
o 

2 A width filter which was considered and rejected for use in our experi- 
ment has the measured characteristic reproduced in Figure 20. In tuning 
o 

to 6840 A, at the blue edge of the Thomson spectrum, the filter shows 
a 70 $ increase in transmitted light for parallel polarization, but a 
40 ^ decrease for perpendicularly polarized light. 

5 .6 System Electronics 

Figure 21 is a functional block diagram of the principal electronic 
components of the system used to fire the laser and the Penning discharge 
and to record Thomson scattering and other data. An asterisk indicates 
those components designed and built especially for this project . 

The laser and the discharge current source are both controlled by 
highly stable delay generators. These allow the time of the laser 
flash to be carefully controlled relative to the plasma evolution and 
also provide control of the discharge pulse duration. Since the dis- 
charge cathodes are directly heated by 60 Hz current of the order of 100 
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Integrated Transmissivity vs Wavelength, Filter 1-2 
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Functional Black Diagram of .System ‘Ele.ct.ronics 












amperes, an alternating axial magnetic field component of peak value 
90 gauss is present in the cathode region, resulting in a significant 
60 Hz modulation of discharge parameters. For this reason, and also 
to eliminate the effect of 60 Hz pickup from other sources, most of 
the experiment is operated from this phase of the a-c line, and all 
experiment events are times to occur within ± 0.05 radians of the zero- 
crossing of this phase. 

The 60 Hz synchronization generator produces a timing pulse whose 
jitter with respect to the 60 Hz zero-crossing is only ± 10 psec or 
± 4 milliradians . By this means, any 60 Hz effect present will certainly 
be quite constant from one shot to the next . 

Of the 3 delay generators, the one provided with the laser intro- 
duces a timing jitter which is at most ± 0 . 05 $ of the delay time or 
±0.5 ps'ec and accuracy of order ± 0.2$ or ± 2 psec. Accuracy of the 
two specially built delay generators is of the same order, but jitter 
is only ± 0 . 005 $ under normal conditions. 

In operation, the synchronization pulse initiates optical pumping 
of the ruby about 1 millisecond before the 60 Hz zero-crossing. When 
the laser lamp begins to conduct, a lamp current rate pulse simultaneously 
starts the 3 timing generators which ultimately start the plasma dis- 
charge current source, open the Pockels cell shutter and then stop the 
discharge, closely concurrent with the zero-crossing. This approach 
allows Thomson measurements to be made at a selected time during the 
plasma discharge without affecting lamp energy input to the ruby rod. 

A vacuum-tube amplifier and high current switch was available for 
the discharge current source. Since the two discharge electron guns 
together are capable of about 600 amperes, temperature limited, it was 
desired to have a current source of similar capabilities in order to 
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get as much plasma electron density as possible. Unfort unatel y, this 
switch is capable of l7o amperes at most; but its other characteristics 
are ideal . It consists of 12 Ediswan I 3 EI beam tetrodes in parallel, 
and when these are driven temperature-limited, they behave like ideal 
current sources. This automatically provides the necessary stabilizing 
factor in the negatiye^resistance region where we operate the discharge. 
In addition, this behavior results in rapid achievement of a steady- 
state discharge in the quiescent mode of operation, and consequent close 
control of discharge input energy, 

Figure 22 is a picture of the typical discharge current and voltage 
profiles, and demonstrates this last point . Other properties of the 
current switch are: maximum open circuit voltage, 9 OOV 5 rise and fall 

times, 5 psec; pulsewidth for 5 $ current amplitude decay, li-00 ps : ec; and 
duty cycle, 0.02 at maximum voltage and current. 

5 .7 Detector Electronics and Photographic Record 

The photomultiplier used to record scattered light signals is an 

RCA 31000F extended red response tube. Twelve stages of electron multi- 

7 

plication provide an overall electrical gain of 6 X 10 at 2350V* Rise 

time of the anode signal is 2.b ns, which is the same as our oscilloscope 

rise time. Although anode dark current is only 50 nanoamperes at this 

voltage and room temperature, the tube is housed in a thermoelectric 

refrigerator. However, the purpose of this device is to eliminate 

temperature-dependent gain changes in the tube rather than to further 

o 

reduce the dark current. We normally operate the tube at 15 C. 

The other major source of gain changes in multiple stage photo- 
multipliers is supply voltage variation since S^nG = njrtn V for an 
n-stage tube. When n ==. 12, we require approximately 0-.1 % voltage 
stability for 1 rfo tube gain stability. Our supply is 5 times better than 


this over a 2b- hour period. 
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Discharge Current 

Vertical: 50 amp/division 
Horizontal: 50 fjsec/division 






Discharge Voltage 
Vertical: 100 volt/division 
Horizontal: 50 fisec/division 



Figure 22 

Typical Penning Discharge Voltage and Current Profiles ( 400 \isec Duration) 
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The photomultiplier " front end" has several special properties which 
are essential for this experiment. The Na/K/Sb/Cs semitransparent cathode 
material has a measured photon-to-elect ron conversion efficiency of 7 - 05 ^ 
at 69^3 A which, at the time of purchase, was the best commercially 
available. This is a factor of 4 improvement over the best standard 
surface (S-20) at this wavelength. High quantum efficiency is extremely 
important from statistical considerations (see chapter 6) when low light 
levels are being detected. Another advantage of the tube is its first- 
dynode gain of 50 , which practically eliminates statistical noise con- 
tributions from the multiplication process itself. Finally, the tube 
is linear up to an anode peak current of I 5 O ma, which corresponds to a 
cathode signal of about 300 photoelectrons with our operating parameters. 
This is 5 times the largest Thomson signal recorded. / 

The photomultiplier anode pulse must be integrated numerically or 
electronically, since it is the total charge which corresponds to the 
scattered photon number. Electronic integration is much the easier of 
the two and also discriminates against high frequency electronic noise. 

A high speed gated integrator was constructed for this purpose and is 
located very close to the PMT chamber, since our laboratory is a highly 
active electrical environment . Its electronic properties are: sensi- 

tivity, 2.4 mv per picocoulomb (pc); time constant, 1 fisec; noise level, 

3 pc; delay, 15 ns; and linearity, about 3 $ U P to Q = 10^ pc or 10 ma 
peak input current, and down to an input pulsewidfh of 5 ns * With our 
operating parameters, the integrator provides an output signal of 22 mv 
per PMT cathode photoelectron. 

Our oscilloscope is a Tektronix r 454. While not the fastest scope 
available at the time of purchase, its 2.4 ns rise time is more than 
adequate for recording the integrated scattering signal, and this input 
bandwidth is maintained down to a sensitivity of 20 mv/division. 
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The oscilloscope camera and CRT together provide a writing speed 
of 3.2 divisions per nanosecond with ASA 10,000 Polaroid film. Features 
contributing to this speed are the lk kv accelerating potential and Pll 
(blue) phosphor used in the CRT, the f^l .2 camera lens and the 0.5 
image-to-object ratio used in this system. 

5 .8 Vacuum System and Mechanical Design 

The vacuum system consists of a mechanical forepump, 300 liter/sec 

oil diffusion pump, liquid nitrogen cold trap and separate experiment 

roughing pump. The cold trap is re-entrant and has a 36 hour liquid 

nitrogen capacity. The vacuum station is quite clean and will achieve 

-8 

a base pressure of 5 X 10 torr. 

The total conductance of the experiment tube is 1 6 liter/sec at 

the vacuum station port, but this is normally reduced to about k liter/sec 

by partially closing the valve to this port, in order to preserve an 

upper limit of 10 $ neutral density variation throughout the discharge. 

The total leak rate of the system is about 0.02 micron-liter/sec, so that 

this pumping speed results in a possible contamination of order 1 $ at 

an operating pressure of 1 p. . Hydrocarbons are responsible for only 15 $ 

of this contamination level, the rest being air or outgassing hydrogen. 

-6 

The tube base pressure is about 2 X 10 torr. Pressures are measured 
at the cold trap port and at the tube fill gas inlet using triode-type 
ionization gages with the appropriate correction factors for different 
gases. These tubes are sensitive to magnetic fields and must be used 
with external field sources turned off. A Granville-Phillips variable 
leak with numerical indicator provides control of gas flow rate into the 
discharge tube. 

The axial magnetic field required for operation of the Penning dis- 
charge is provided by 8 air-core water-cooled coils made by Ogallala 
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Electronic Mfg Inc. Each has a 20.3 cm bore, overall diameter of 38 . 1 
cm and 8.0 cm thickness, and is rated to produce a maximum field of 
1200 gauss on axis, at 16.5 gauss/ampere . With a 15*3 cm center-to-center 
spacing, the assembly produces about 30 gauss per ampere on axis with 
a maximum spacial variation of ±2 $. The available current supplies for 
these coils limit the available field intensity to 750 gauss, but more 
is not needed for our purposes. 

In the laser beam optics, angular alignment of the order of 5 
seconds of arc must be maintained to assure proper registration of the 
beam focus with the observation volume . This is just the price we pay 
for using distance to reduce equipment scattering. Prevension of 
unintended motion of the scattering volume within the plasma column 
depends on rigid interconnection of the Penning discharge and the experi- 
ment optics. A nonmagnetic bench was constructed for this purpose from 
structural aluminum. The main beams of the bench are sufficiently 
rigid to deflect less than 0.001 inches per 800 pounds of distributed 
load . All other elements of the bench which might affect optical align- 
ment, including lens and mirror holders, are also aluminum. In this way, 
misalignment due to differential thermal expansion is avoided. 

Operational alignment of the entire system is provided by a small 
gas laser beam^ expanded to the same size as the ruby beam and arranged 
so that insertion of a right angle prism produces a test beam which is 
concentric and collinear with the ruby laser output. This beam was also 
used for initial alignment of the detection optics, by using a small 
mirror to redirect it from the center of the scattering volume. 
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6. RESULTS AND EXPERIMENTAL PROCEDURE 


In this chapter, we present and interpret our Thomson scattering 
measurements, and discuss relevant data obtained by other means. Experi- 
mental conditions are analyzed carefully and related to the procedures 
used to collect and analyze data. 

6 ,1 Langmuir Probe Measurements 

When Thomson scattering measurements are possible, they give much 
more detailed and reliable information about the state of a plasma than 
do Langmuir probe measurements. However, the probe was invaluable as 
an exploratory device for locating the best combination of ioniza- 
tion fraction, electron density and parameter repeatability from shot 
to shot . Even if these measurements had given different temperatures 
from those deduced from Thomson measurements, we would not have been 
too upset, because probe electron collection is strongly influenced 
by magnetic fields. However, it is gratifying that good agreement was 
found . 

To take Langmuir data, the time-sample technique introduced by 

86 

Hosea was used. This method is desirable in a pulsed plasma where 
density and temperature change significantly during the plasma evolution, 
and permits rapid construction of a probe curve characteristic of a 
particular time during the pulse. In this technique, an oscilloscope 
is set up to record probe voltage versus current continuously, but 
is gated "on tf for only a brief interval (0.5 psec) during each plasma 
pulse. With a stable plasma firing repeatedly, the result is a fixed 
blinking dot on the scope screen, which will trace out a probe curve 
representing the selected time epoch when probe voltage is varied. 
Further, the extent to which the blinking dot remains confined to a 


82 




single curve in the V-I plane is a good quantitative measure of the 
reproducibility of plasma density and temperature. 

Typical probe traces taken by time sampling are reproduced in 
Figure 23- These pictures represent conditions during the discharge, 
for magnetic fields of 150 and 600 gauss, respectively. The transi- 
tion to an unstable discharge is clearly seen. In the second picture, 
neither the space nor the floating potentials are defined, and even 
an approximate temperature cannot be deduced. However, Thomson measure- 
ments have been made in a similar environment (series 4). 

A summary of the variation of plasma density and temperature with 
changing neutral pressure at a time 25 ^sec before shutdown of the dis- 
charge current source is shown in Figure 24, as measured by the Langmuir 
probe. These measurements were made at the center of the plasma column 
with a uniform magnetic field of 190 gauss. The subsequent variation 
of density and temperature indicated by the Langmuir data for the same 
magnetic field and radial position was shown in Figure 5* Curves for 
other pressures in the range shown in Figure 24 behave in very much 
the same way in the afterglow. Below this pressure range, the available 
neutral number begins 'to limit the electron density while above it, the 
ionization fraction is so small that the electron temperature falls 
below useful values for our experiment. 

6.2 Rayleigh Scattering 

The total Rayleigh scattering cross-section is 


dCT R 

an 




sin 


2* 


(6-1) 
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B = I 9 O gauss 
(Quiescent) 




B = 600 gauss 
(Unstable) 


Figure 23 

Langmuir Probe Curves Showing Effect of Strong Instabilities 
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in a gas of neutral density n^ and refractive index ja, for radiation of 
wavelength \ Q • Doppler shifting is present here, too, of course, but 

the molecular velocities are so small that the thermal doppler half-width 

° 3 o 

can be no more than, say, 0.1 A for a molecular temperature of 10 K. 

Accordingly, in a Thomson scattering environment, the Rayleigh scattering 

can almost always be thought of as occurring at the laser frequency only. 

The Rayleigh cross-section in the scattering plane (sinY = l) is 

“28 2 —28 2 
2.071 X 10 cm for nitrogen and only 0.4^9 X 10 cm for hydrogen, 

at 69^3 A. The ratio do /do of the total Thomson to the total Rayleigh 

T R 

cross-sections is therefore 1,720 per atom and electron in hydrogen. The 
relative magnitude of the observed Thomson and Rayleigh signals at the 
laser line will depend on the ionization ratio and on the fraction of 
the total Thomson spectrum admitted by the dispersing instrument. A 
more useful figure is our measured Rayleigh scattering contribution of 
about 0.4 photoelectrons per joule-micron (light energy, neutral pressure). 
In a typical Thomson measurement, our ionization fraction was such that 
the Rayleigh contribution was no more than 10$ of the equipment- 
scattered light at the laser frequency, or 0.5$ tlie Thomson signal 
at the spectrum center. 

Therefore, Rayleigh scattering is not an important consideration 
in analyzing our Thomson results, but it is an extremely useful way of 
calibrating and fine-aligning the light detection system, as well as 
checking for drifts in the original filter calibration. For example, 
through Rayleigh calibration, we find no more than a factor of 2 total 
difference between calculated and measured absolute scattering magnitudes, 
a figure which includes absolute error in pressure indicated by our 
ionization gages. 

Rayleigh measurements were made at a pressure of 500 m. in d r Y 
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nitrogen in order to get a signal of the order of 10 photoelectrons 
at pressures low enough to avoid photoionization of the gas. Pressure 
was established by a ratio method, at flow rates of order 0.2 liter/sec 
in the discharge tube. 

6. 3 Gas Selection 

Hydrogen was used for all Thomson measurements to be reported. 
Among the gases which could be used in a Penning discharge, Table 6-1 
lists the best nine in decreasing order of desirability, with regard 
only to their spectral activity in the Thomson region. 

Table 6-1 

Spectral Activity of Various Gases in the 
Interval 6700 $ - 7100 8 , from M.I.T. 

Wavelength Tables 


gas 

number of 
spectral lines, 

, 6700-7100 & 

location of 
nearest lines 

reactive? 

easy to use? 

1 . hydrogen 

0 

+1,261, -380 

no 

yes 

2. sodium 

0 

+ 1 , 240 , -782 

yes 

no 

3 . helium 

2 

+ 123 , -265 

no 

no 

4. rubidium 

3 

+99; -138 

yes 

no : 

5 . nitrogen 

17 

+3, -0 

no 

no 

6 . oxygen 

5 

+ 59 ; -32 

yes 

no 

7 . chlorine 

11 

+ 8 , ^9 

yes 

no 

Q'. potassium 

4 

+21, -5 

yes 

no 

9 . neon 

8 

+81, • -13 

no 

yes 


The spectrum of argon contains 35 lines in this interval, although 
it is a very well-behaved gas for Penning discharges. Either hydrogen 
or helium might have been used, except that helium is more difficult 


to use in the discharge. 
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6 . k Plasma Heating by the Laser 


The equivalent fractional increase in electron temperature which 
can occur via inverse bremsst rahlung during the laser pulse is given by 

bl 

Kunze as: 


-7 n i Z g ff X c 


- e _hv o /kT e 


where is electron temperature expressed in ev , n^ is ion 

density, J is laser energy intensity in the scattering region (joules/ 
cm ), is the laser wavelength in cm, and g^^ is a Gaunt factor 

which may be approximated by unity for purposes of rough estimation. 

69 

This expression is just Spitzer' s formula for the inverse 
bremsstrahlung absorption coefficient for a fully ionized plasma. How- 
ever, it includes the assumption that electron-electron relaxation times 

t in the plasma are several times smaller than the laser flash dura- 
ee 

tion t , which is in turn supposed to be several times smaller than 
L 

the energy equilibration time t^ for plasma electrons and ions as 

species. That t < t in our plasma can be demonstrated with our 
ee L 

handy formula for this quantity, also based on Spitzer: 


q/2 

t ^ 300 /n _ ^ nanoseconds . 
ee ^ e e !2 


(6-3) 


12 

If, for example, = 3 ev and n e \2 = ^ ( 10 cm J ) , typical values, 

we find t ~ 7-8 nanoseconds, so that the usual laser flash persists 
ee 

for 4 or 5 electron relaxation times. Spitzer' s formula for t gives 
a time of several microseconds, but the actual time will be closer to 
hundreds of microseconds in our discharge because the hot ion gyroradii 
are so large that they instantly escape. 
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With these provisions, for at hv = 1.79 ev > Kunze's 

formula can be simplified to 


A0 e -7 

T ,u,w 

e 



(6-M 


2 

which with our typical energy intensity J = 150 joules/cm , ©^ = 3 ev 

and n el 2 = ^ gives an estimate of the relative heating A® 0 /® e ~ 

-5 50 

5 X 10 . It should be noted that the approach used by Gerry^ lead's 

to an unnecessarily large estimate for the relative heating. 


6.5 Plasma Light 

For our hydrogen plasma, bremsstrahlung is expected to be the most 
important process producing plasma light in the Thomson interval. We 
employ a formula of Glasstone and Lovberg for the bremsstrahlung produc- 
tion rate, converted to practical units for a hydrogen plasma and for 



V x 


* 3-9 x 10 


■12 


2 - 
n el2 S 


- 1 . 79 /©, 


, O “5 

e dA watts/cm J 


(6-5) 


into sterradians. Since our Thomson data will be expressed in units 
of PMT cathode photoelectrons and since we have not yet taken into 
account- the light collection efficiency of our detection system, a more 
useful form of (6-5) for direct comparison with measurements is, again 
writing g = 1: 


2 -1/2 -1 7 Q /0 0 

I as 1.1 X 10 0 n, © e ’ {y/ * Tl] dQ dV dA 

p el2 e 

photoelectrons/nanosec . 


( 6 - 6 ) 


89 



In (6-6), T is the effective filter transmissivity, equal to its abso- 

o 

lute integrated transmissivity divided by its bandwidth dA in 
Angstroms, dQ is the viewing solid angle, ^ is the PMT quantum effi- 
ciency and dV is the visible plasma volume (larger than the scattering 

o 

volume) in cm . Inserting the correct values for these quantities, we 

find for n _ ^ = 20 and © = 3 ev that I = 3 X 10 ^ p.e./ns. This 

el2 e P 

calculation agrees well with our measurements, which show less than 

-3 

2 X 10 p.e./ns from all sources including the PMT dark current noise 
discussed in the next section. 


6 .6 Photomultiplier Produced Noise 

The recorded variable is the quantity of charge contained in the 
PMT anode pulse. Because the PMT has unusually high first dynode gain, 
the largest single source of statistical or "shot" noise is in the con- 
version of photons to photoelectrons at the cathode, rather than in 
electron multiplication. 

Accordingly, we expect that since the arrival of N photons consti- 

P 

tutes N Bernoulli trials in which "success" is the production of a 
P , 

photoelectron with probability 0 . 07 j the photoelectron statistics as 
measured at the anode will be governed by the binomial distribution, with 

0.07: 

( N \ N (N -N ) 

N p |n e (i-Ti) p e . (6 -t) 


Further, we expect that the normal density will be a good approximation 

to the binomial distribution for most of our conditions, so that the 

expected standard deviation of the number of photoelectrons received 

due to the arrival of N photons will be 

P 

¥ 

This is the photocathode quantum efficiency. 
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y 


(6-8) 


a = IWh) « \Tnti = 'liT 

where is the total number of photoelectrons received, whether in 

a single shot or over a number of shots. Accordingly, at regions of 
very small signal, a sufficient number of shots is taken to reduce the 
standard deviation to a reasonable level. 

With no light input, a photomultiplier cathode or dynode produces 
a very small thermal current which is called the M dark current” when it 

is measured in amplified form at the PMT anode. At the tube operating 

o , 

temperature of 15 C, and ah operating voltage of 2359 V, the PMT dark 

current is only 50 nanoamperes. Referred back to the cathode, the equi- 

“6 

valent dark photoelectron production rate is then only 6 X 10 p.e ; /ns, 
which is entirely ignorable for our purposes. 

6.7 Plotting Thomson Data 

We are very fortunate in having designed an experiment in which the 
statistical noise expressed by (6-8) is the only significant source of 
noise which must be dealt with over the scattered spectrum. Of course; 
equipment-scattered light including Rayleigh scattering is present at 
the spectrum center, and must be subtracted. In general, this has been 
measured repeatedly in each Thomson series, and its standard deviation is 
very low compared to that of the signal. In any case, the amount sub- 
tracted is plotted along with each scattered spectrum reported, and 
labeled "ESL" . 

To take data, Langmuir measurements were initially employed as des- 
cribed earlier to locate an operating point giving sufficient density 
and temperature to allow measurements. The plasma discharge was then 
fired recurrently at about 1 pps for 3 or 4-, hours to permit sputtered 
impurities to be flushed out of the system. This could be verified when 
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no change in tube pressure resulted from firing the discharge. Then, 
Thomson measurements were made at each wavelength listed on a precalcu- 
lated shot schedule. This schedule was arranged so that approximately 
7 measurements would be made on each side of the spectrum, including at 
least 2 at the center and 2 in the extreme wings, for each of the two 
scattered spectra (\\ and ±) to be measured in one series. The schedule 
sequence begins at the spectrum center and moves through the blue wing, 
then returns to the spectrum center and moves through the red wing, always 
alternating between a point for the parallel and a point for the perpen- 
dicular spectra. Each "point" in these measurements consisted of from 
10 to 20 repeated shots, the number required depending on received signal 
strength and our estimate of the desirable standard deviation. 

Since this procedure already sounds like the tale of the man going 
to St. Ives, we might mention that a typical M series”, leading tq one curve 
each for the parallel and perpendicular spectra, required approximately 
1,000 shots, 40 rolls pf Polaroid film and 10 hours to complete. This, 
in fact, is the reason for the complicated procedure listed above. We 
trust our measurements because the same signal was found in all cases 
upon returning to the spectrum center after 5 hours spent in the blue 
wing. Some examples of scattering data are reproduced in Figure 25 • 

Later, the signals were read off the Polaroid prints, numerically 
averaged for each measurement wavelength, and corrected for ESL and for 
changing filter integrated transmissivity and laser output energy. 

The resulting values, normalized to normal incidence transmissivity and 
1 joule of light, are plotted on a linear scale versus [ AX/sin( 9/2) ] . 

After finding the wavelength around which the plotted values were symme- 
trically distributed, the logarithm of signal strength was replotted 
versus .[ A\/sin( 9/2) ] . It will be seen from equation (2-30) that such 
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Scattered Energy Signal 
(2 shots) 

Vertical: ~ 5 photoelectrons 
per division 

Horizontal: 20 ns/division 



Laser Energy Monitor 
Signal (5 shots) 

Vertical: 0.2 joules/division 
Horizontal: 20 ns/division 


Figure 25 

Examples of Scattering Data 
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if the 


a curve will lie on a straight line of slope equal to ~(a® e ) 
electrons are Maxwellian. For a plot on decilog paper, a = 1738. 

In fact, all our measurements show a Maxwellian plasma, except the 
last, for which the discharge was deliberately operated above the cri- 
tical field for instability. Even for this case, the data conversion 
method used here was the most appropriate. 
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7, ACTUAL THOMSON SCATTERING RESULTS 


7.1 Data Presentation 

A total of four measurement series were completed, each consisting 
of a parallel and a perpendicular distribution measurement. All were 
made at the plasma column center. The first series, however, featured 
a smaller sampling volume than the last three. For series 1, the sampling 
cylinder is 3*0 mm and 5.3 mm in length for the longitudinal and trans- 
verse parts, respectively. This is the only complete series in which 
the x and || data were not interleaved. Series 2 , 3 and 4 were made with 
a cylinder length of 8.5 mm. The sampling cylinder has a diameter of 
1 mm, determined by the laser beam waist. Total plasma diameter along 
the cylinder direction is more than 4 times the cylinder length. 

Discharge time for all measurements is 25 1 -isec before the cessation 
of discharge current. Magnetic field strength is a uniform 190 gauss 
for series 1-3 and 450 gauss for series 4, where we deliberately ex- 
ceeded the critical field for plasma instability. Neutral pressure varies 
among the series, and is listed for each set of figures. 

Series 1 was our very first successful Thomson measurement (Figures 
26-28). At that time, we thought we would have to be satisfied with a 
signal strength of 5 -IO photoelectrons per shot, and so completed our 
only high temperature measurement. In succeeding measurements, this was 
never attempted, since high temperatures go with low densities and the 
two effects together reduce signal strength to masochistic levels. While 
the standard deviations achieved in this series are not impressive, 
better results would most likely have been impossible with our apparatus. 
Within experimental error, series 1 shows two equal temperature Max- 
wellians with a temperature of about 5 ev. The most extreme point on 
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Series 1, Azimuthal Distribution 
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Series 1, Longitudinal Distribution 
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Figure 29 

Series 2, Longitudinal Distribution 
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Series 2, Logarithmic Plot 


the transverse plot is significant because of its very low magnitude, 

0.15 db 0.15 p.e./shot. This point represents an average incidence 
of 2 photons per shot. 

Series 2 is an isolated longitudinal plot (see Figures 29, 30) which 
could not be accompanied by a transverse measurement due to laser failure. 
It is included because it is the first instance of what we consider good 
standard deviation in the results. The accompanying log plot (Figure 30) 
shows that a single temperature of 3.18 ev can describe all points on 
this curve, out to electron energies of about 13 ev . 

Series 3 (refer to Figures 3 ^ "33) our bes * overall pair of measure - 
ments, showing excellent standard deviation and little if any serious 
evidence of the high energy tails to be ■ found in series 4. Single tempera - 
tures describe all points on both curves out to electron energies as 
large as 6 times the thermal energy, and these temperatures are equal 
within experimental error. These temperatures are 3.0 1 0.3 ev for the 
perpendicular distribution and 3*2 ± 0.2 ev for the parallel distribu - 
tion. The measurements on these plots show Maxwellian behavior over a 
much broader range than has ever been seen before in Thomson scattering . 

Series 4 (Figures 3^“37) provides by far the most interesting of 
our results. It is immediately clear from either log plot that a straight 
line will not fit the data within experimental error. In fact, the dis- 
tributions have quite large high energy tails, so that a free-form curve 
drawn through the data points has a tangent line at zero energy corres- 
ponding to about 2 ev, while in the wings of the distributions, the 
tangent line slope is 15 to 20 ev in both cases. It is worth noting 
that the transverse plot has lower " central temperature" and higher "wing 
temperature" than the longitudinal plot. 

The accompaning gaussian curves on the linear plots (Figures 34 
and 36 ) for the two directions correspond to a temperature of 2.27 and 
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Figure 32 

Series 3, Longitudinal Distribution 
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Series 4 , Longitudinal Distribution 
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Series k, Logarithmic Plot of Longitudinal Data 
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Series 4, Azimuthal Distribution 
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Figure 37 

Series 4, Logarithmic Plot of Azimuthal Data 



2.32 ev, respectively. These were chosen as the best fit to the low 
energy data points in each log plot, and are meant to demonstrate the 
excessive electron numbers in the wings. However, the log plot is 
more dramatic for this purpose. 

This is the first time that measurements showing this behavior have 
been made using Thomson scattering. Of course, the standard deviations 
are much worse in series k because of the chaotic conditions in the 
plasma they represent. Nevertheless, one should keep in mind the near 
impossibility of making measurements in this regime by any other means . 

To put this point in perspective, we again refer to Figure 22, which 
shows a Langmuir trace taken under the same conditions . 

7.2 Data Summary and Suggestions for Future Work 

These measurements demonstrate nearly identical behavior in the 
electron velocity distributions parallel and perpendicular to the magnetic 
field direction in a highly ionized Penning discharge in hydrogen. Even 
when the discharge is operated so as to encourage the Hoh instability, 
the same general behavior is found between the distributions, with the 
exception that the azimuthal distribution seems to have the lower central 
temperature and the higher wing temperature of the two. Otherwise, when 
a Maxwellian exists (series 1,2 and 3 ), the temperatures found are equal 
to within experimental error and agree well with Langmuir probe measure- 
ments taken under the same conditions. 

These results are quite different from those found by Koons and 
Fiocco^'^ but, as explained earlier, the conditions of their experi- 
ment were quite different from our own. On the other hand, our results 

6b 

agree qualitatively with those being obtained by McCormick ^ on a highly 
ionized hollow cathode arc discharge in argon. 
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We are not surprised by the findings. On the contrary, the electrons 
should thermalize very rapidly in this discharge at points relatively 
far from the chaotic cathode fall region near the electrodes. 

Qualitatively, this experiment is the best yet reported in the 
following areas: 

(1) Equipment scattered light is several orders of magnitude below 
other comparable experiments, due to careful design of stray light 
trapping systems. 

(2) This is the first Thomson experiment to measure Maxwellian 
behavior over a range from zero to 6 thermal energies. 

(3) This is the first Thomson experiment to clearly show non- 
Maxwellian behavior in a Penning discharge plasma. 

Also, it can be deduced from the data that good k-space resolution 
can be had in a Thomson scattering measurement, and yet permit measure- 
ments out to at least one thermal energy for electron densities as low 

12 -3 11 -3 
as 10 cm . Of course, densities down to 10 cm or lower could be 

permitted in an isotropic plasma where light collection from many direc- 
tions is possible. The idea that Thomson scattering may not be used 
13 -3 

below 10 cm is just a myth. 

Up to this point, we have not commented on the apparent slight zero 
offset toward the blue shown in most of our data. It is not clear to us 
if it is significantly larger than the experimental error, since this 
error could arise from so many different parameters. Nevertheless, where 
it occurs, it always occurs in the same direction, and amounts to only 

o 

3 or 4 A. It could be due to net electron drift, since it occurs with 

one exception (series 1, azimuthal) in the longitudinal distribution. 

If this is the case, the magnitude of such a drift would be of order 
7 

10 cm/sec . 
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We note here that this technique is completely capable of demon- 
strating plateaus or "bump in tail" effects in the electron velocity 

88 

distribution, such as discussed in the theory . It would be particu- 
larly interesting to extend Thomson measurements to much lower densities, 
in order to observe longer-lived behavior of this sort than we have 
found. Such experiments could be done, since at the very low densities 
the plasma heating limit is much less severe (see 6-4), allowing the 
use of very high laser energies provided that ( 6 - 3 ) is observed and 
the laser power kept constant. For example, a plasma with temperature 
3 ev but density 10*^ cm ^ could be observed by our optics using a 10^ 
joule laser with a 30 microsecond pulse length, without altering the 

heating fraction of 5 X 10 which we obtain. Lower densities yet could 

2 

be observed by increasing the light collection solid angle by 10 to 1 

3 

sr. One must be careful, however, to maintain n D > 1 unless consi- 

e e 

deration is given to the consequences. 
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